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■ Cosmological perturbations are generally described by quantum fields on (curved 

' but) classical space-times. While this strategy has a large domain of validity, it can 

CN ■ not be justified in the quantum gravity era where curvature and matter densities are 

C , of Planck scale. Using techniques from loop quantum gravity, the standard theory of 

\ cosmological perturbations is extended to overcome this limitation. The new frame- 

\^ I work sharpens conceptual issues by distinguishing between the true and apparent 

trans-Planckian difficulties and provides sufficient conditions under which the true 
difficulties can be overcome within a quantum gravity theory. In a companion pa- 
per, this framework is applied to the standard inflationary model, with interesting 
implications to theory as well as observations. 



I. INTRODUCTION 



(N 
> 

^ ■ This is the second in a series of papers whose goal is to investigate whether current sce- 

narios of the early universe admit quantum gravity completions and, if they do, to study 
the implications of the resulting Planck scale dynamics. The first paper [1] summarized 
CN ' the underlying framework and its applications to inflation for a broad audience of theo- 
retical physicists. This paper is primarily addressed to the quantum gravity community 
and provides a detailed extension of the cosmological perturbation theory to the Planck 
^ \ regime. Specifically, we will consider gravity coupled to a scalar field and study the dynam- 
^ ■ ics of quantum fields representing scalar and tensor perturbations on quantum cosmological 

■ ■ ■ space-times. In the third paper [2], addressed to cosmologists, this framework is used to 

show that the inflationary scenario admits a quantum gravity extension and to analyze the 
physical implications of pre-inflationary dynamics, i.e. the quantum evolution from the big 
bounce of loop quantum cosmology (LQC) to the onset of the standard slow roll inflation. In 
the future, we hope to examine whether alternatives to inflation also admit viable quantum 
gravity completions and, if so, explore the resulting Planck scale physics. 

In the theoretical explorations of the early universe, one generally uses the Friedmann, 
Lemaitre, Robertson, Walker (FLRW) solutions to the Einstein equations (with appropriate 
matter sources) as background space-times. The focus is on the dynamics of quantized fields 
representing linearized perturbations propagating on these backgrounds. (See, e.g., [3-6].) 
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The necessary framework of quantum theory of hnear fields in curved space-times has been 
well developed, thanks to the ongoing research that date back to the mid 1960s. (See, e.g., 
[7-9]). However, the FLRW space-times of interest are invariably incomplete in the past 
due to the big bang singularity where matter fields and space-time curvature diverge. It is 
widely believed that general relativity is simply not applicable once curvature reaches the 
Planck scale, whence there is no justification for using quantum field theory on solutions 
to Einstein's equations in this domain. Quantum gravity must intervene in an important 
fashion. Thus, to encompass the Planck regime, one needs a quantum gravity extension of 
the standard cosmological perturbation theory. 

Loop quantum gravity (LQG) provides a promising avenue to meet this goal because by 
now the big bang singularity has been resolved in a variety of models in LQC; the k=0 FLRW 
model on which we will focus [10-15], as well as their generalizations that include spatial 
curvature [16, 17], a cosmological constant [18-20], anisotropies [21-24] and the simplest 
type of inhomogeneities and gravitational waves [25-29]. (See, e.g., [30, 31] for summaries 
of these developments.) It is therefore natural to use LQC as the point of departure for 
extending the cosmological perturbation theory. However, to do so, we cannot just mimic 
the standard procedure used in general relativity because LQG does yet offer the quantum 
version of full Einstein's equations which one can linearize around a quantum FLRW space- 
time. Therefore we will use the general strategy that has been repeatedly followed in LQG: 
First truncate the classical theory in a manner appropriate to the physical problem under 
consideration, then carry out quantization using LQG techniques, i.e., paying due attention 
to the underlying quantum geometry, and finally work out the physical implications of the 
framework. This strategy has led to advances in the cosmological models referred to above, 
as well as in the investigation of quantum black holes [32-34] and the spin foam derivation 
of the graviton propagator [35-37]. 

To extend the cosmological perturbation theory, then, we will begin by focusing on the 
following sector of the full phase space of general relativity: homogeneous, isotropic con- 
figurations together with first order inhomogeneous perturbations. However, to encompass 
the Planck regime, we must now use a quantum FLRW geometry as background, and study 
the dynamics of quantum fields representing scalar and tensor modes propagating on these 
quantum geometries. Since such a quantum geometry provides only probability amplitudes 
for various FLRW metrics, we no longer have a sharply defined, proper or conformal time. 
How can one then describe the dynamics of inhomogeneous perturbations? In [38] , this issue 
was resolved for test quantum fields on quantum FLRW space-times by deparameterizing 
the Hamiltonian constraint in the background, homogeneous sector. Then one can regard 
the background scalar field as a relational time variable with respect to which physical 
observables evolve. This is a new conceptual element, made necessary by quantum gravity 
considerations. We will use the same strategy. 

However, to encompass cosmological perturbations, we will need three significant exten- 
sions of Ref. [38]. First, while that work discussed a test quantum scalar field, now the 
test fields include metric perturbations. Second, to systematically arrive at the evolution 
equation of perturbations on quantum geometry, one needs an improved strategy. The third 
and most important difference is that, because of its focus on conceptual issues such as the 
problem of time, the analysis in [38] was restricted only to a finite number of modes of the 
test scalar field, and thus avoided the ultraviolet difficulties from the start. In this paper, 
by contrast, we do not truncate the number of modes and much of the difficult analysis is 
devoted to these ultra-violet issues. 
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Indeed, these issues play a central role in testing self-consistency of our procedure. The 
key approximation underlying our truncation strategy is that inhomogeneities can be re- 
garded as perturbations — i.e., their back-reaction on the space-time geometry can be ig- 
nored. In the classical theory, a solution obtained using this approximation is regarded as 
self-consistent if the stress-energy in its inhomogeneities is negligible compared to that in 
the homogeneous background for the entire dynamical regime of interest. In our quantum 
theory of the truncated phase space, the situation will turn out be the following. Fix a quan- 
tum FLRW background geometry, described by a state \l/o, which evolves (with respect to 
the deparameterized, internal time) via a Hamiltonian Ho- The state undergoes a quantum 
bounce at some time = 0^. There are natural initial conditions for quantum states ip of 
perturbations at = 0g. We will show that for these states, it suffices to focus just on the 
energy density, rather than the full stress-energy. If the inhomogeneities are to be regarded 
as perturbations, we have to choose the initial if) so that the energy density in it is small 
at = </)g. The question is whether it continues to remain small in the entire duration of 
evolution of interest. This is not at all guaranteed, especially because of the Planck scale 
curvature during and following the superinflation phase immediately after the bounce. But 
if it does, then '^o®'4^ would be a self-consistent solution to the truncated quantum theory.^ 
Note that the very formulation and subsequent analysis of this self-consistency criterion 
requires a well-controlled definition of the Hamiltonians Ho and Hi (and the corresponding 
energy density operators). Existing results on LQC already provide a well-defined, specific 
Ho- On the other hand. Hi is a composite operator on the Hilbert space of perturbations; its 
formal expression involves products of operator valued distributions. Therefore, an appropri- 
ate ultraviolet regularization and renormalization is essential, first to obtain a well-defined 
evolution of ip, and second to check that the energy density in the state if) continues to 
remain small compared to that in the background. To properly handle this issue, we will 
carry over the well developed techniques of adiabatic regularization from the quantum field 
theory on classical FLRW space-times to that on quantum FLRW space-times. 

Because this article is addressed primarily to the quantum gravity audience, we will make 
two assumptions that will enable us to have good mathematical control without unduly 
simplifying the essential conceptual underpinning. First, we will assume that the spatial 
topology is that of a 3-torus rather than M?. (In practice, there would be no obvious 
conflicts with CMB observations if one lets the physical radius of each of the three S*^ in T'^ at 
the last scattering surface to be greater than the known radius -Robs of the observable universe 
at that time [39].) This assumption lets us cleanly avoid spurious infinities that would arise 
with an M? topology simply because the background fields are homogeneous. However, our 
main results extend to the topology and at various junctures we will indicate how these 
issues are handled in that case. Second, we will assume that there is no potential; V{(f)) = 0. 
To incorporate inflation, one has to remove this assumption and the necessary modifications 
are summarized in the Appendix A of Ref. [2] which discusses inflation in detail. 



Of course, sclf-consistcncy by itself docs not imply that a truncated solution is necessarily close to an exact 
one. This is so even in classical general relativity because, even if the first order perturbation remains very 
small, the sum of all higher order terms could be large. Yet, in classical cosmology and investigations of 
black hole perturbations, if the back-reaction due to first order perturbations remains negligible — i.e., if 
the test field approximation is self consistent — the first order truncation is generally regarded as a good 
approximation. We adopt the same viewpoint in this paper. 
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The paper is organized as follows. To put this work in a broader perspective, in section II 
we first summarize a few procedures that have been followed in the literature to extend the 
cosmological perturbation theory using LQC, point out their merits and limitations, and then 
present our strategy. In section III we spell out the desired truncation in the classical theory. 
Specifically, we start with the full phase space of general relativity coupled to a scalar field, 
truncate the constraints to second order in inhomogeneities, and express the Hamiltonian 
Hi governing dynamics of the first order perturbations in terms of those gauge invariant 
variables that are most convenient for passage to the quantum theory. In this passage, 
certain conceptual subtleties arise that affect subsequent technical results. We discuss them 
in some detail because, while they are well known in the discussion of perturbations of black 
holes, they are often overlooked in the cosmological context, especially in the LQG literature. 
(See Appendix A for a discussion of these issues in a simpler context of the theory in 
Minkowski space-time.) 

In sections IV - VI we discuss the quantum theory. The main steps involved in the 
construction are first collected in section IV; the detailed construction of the Hilbert space 
of states follows in V; and the necessary regularization and renormalization of composite 
operators is then discussed in section VI. 

Physical states 0; Q, T, ) depend on i/ ~ (the cube of the scale factor), the back- 
ground scalar field 0, and on the gauge invariant scalar mode Q and the two tensor modes 
T of inhomogeneous perturbations. As is standard in LQC, the form of (the homogeneous 
part of the) quantum constraint will enable us to regard as a relational time variable so 
that, at a fundamental level, dynamics refers to this internal or emergent time. Since the 
back-reaction of perturbations is neglected, one can express \E'(z/, 0; Q, T) as a tensor product 
\l/o(z/, 0) ® -0(2, T, 0) and study its evolution. Then, a surprising and key simplification oc- 
curs: The evolution of the quantum fields Q, T on the quantum geometry defined by \E'o (in 
relational time 0) is mathematically equivalent to their evolution on an effective background 
metric cjab (in its conformal time fj) where Qah is an effective metric, 'dressed' with quantum 
corrections. Physically, while the evolution of the quantum fields Q, T is sensitive to the 
quantum nature of geometry defined by it does not 'see' all the details of this quantum 
geometry: It is sensitive only to the expectation value and certain aspects of the fluctuations 
of the quantum metric which, it turns out, can be captured by gab- At late times (jab can 
be well approximated by a FLRW solution to Einstein's equation. Therefore, this passage 
from the relational time to f/ also makes it manifest that, while the quantum constraint 
provides an evolution w.r.t. starting right from the bounce, away from the Planck regime 
this evolution reduces to that used in the familiar treatment [3-6] of quantum perturbations. 

The key question, as we have already indicated, is whether our basic assumption that 
the back-reaction can be neglected is borne out in the resulting solutions, especially in the 
Planck era. The technical discussion of section VI is devoted to this issue and leads to a 
sufficient condition for the solution \E'o(i^, 0) ®'0(Q, T, 0) to be self-consistent within the 
truncation approximation. In the next paper [2] we will show that for a large class of initial 
conditions, this criterion is met in presence of a quadratic potential V^(0) all the way from 
the bounce to the onset of slow roll inflation. In section VII we summarize the main results 
and discuss various issues, including the choice of initial conditions. 

Our conventions are as follows. The space-time metric will be assumed to have signature 
-,-!-,+,+. we set c = 1 but keep {G and) h explicitly in various equations to facilitate 
the distinction between classical and quantum effects. Finally, as in the quantum gravity 
literature, we will use Planck (rather than the reduced Planck) units. (Thus, our £pi = yJhG 



5 



and our Planck mass mpi = ^/h/G is related to the reduced planck mass Mp\ via mp\ = 

II. STRATEGIES INSPIRED BY LQC 

The singularity resolution in LQC has motivated a large number of investigations aimed 
at incorporating the underlying quantum geometry effects in the standard cosmological 
paradigms. In this section we will first briefly summarize the main strategies and then 
explain the approach that is followed in the rest of the paper. This concise summary should 
help non-experts to see the inter-relation (or lack thereof) between the rich set of ideas that 
are being pursued in the LQG literature. For experts, it should help clarify some subtle 
issues that have not been emphasized in the literature and also bring out the continuity and 
coherence of the approach used here. 



A. A brief summary 

Recall first that in the FLRW models the LQC quantum states of interest remain sharply 
peaked along the bouncing solution of effective equations that incorporate leading order 
quantum corrections. This surprising behavior was first encountered in numerical simula- 
tions but subsequent work led to an analytic understanding through a change of represen- 
tation [14], use of coherent states [40] and the WKB approximation [41]. Early attempts 
sought to exploit this property to incorporate quantum gravity corrections right before the 
onset of the slow roll inflation [42], or, for the evolution of perturbations during inflation 
(see, e.g., [43]). However, to incorporate these corrections to the evolution, the LQC effec- 
tive equations for the FLRW background do not suffice; one also needs the LQG modified 
perturbation equations. Since a well established set of quantum Einstein's equations is not 
yet available in full LQG, the strategy was to simply use the first order perturbation equa- 
tions from general relativity, the FLRW background space-time being simply replaced by 
the effective solution in LQC. Unfortunately, this procedure is conceptually unsatisfactory 
because the background LQC space-time is no longer a solution to Einstein's equations. 
Indeed, now there is an ambiguity in what one means by 'linearized Einstein's equations': 
Two sets which are equivalent when the background is a solution to the exact Einstein's 
equations are generically no longer equivalent if the background is a solution to a different 
set of equations. To our knowledge, a systematic procedure to handle this ambiguity was 
not part of the general strategy. 

Another set of papers, geared to capture the 'inverse triad (or inverse volume) corrections' 
of LQC, is based on the notion of 'lattice refinement' (see, e.g., [44-46]). In the homogeneous, 
spatially compact case, these corrections are meaningful and physically interesting. In the 
spatially non-compact case, typically the corrections depend on the choice of a fiducial cell 
used as an infrared regulator [13] and disappear when the regulator is removed [31]. In more 
recent versions, one considers spatial topology but decomposes the spatial manifold in 
elementary cells and approximates the inhomogeneous configurations of physical interest by 
configurations which are homogeneous within any one cell but vary from one cell to another. 
Physically, this is an attractive strategy. However it requires a fresh input — that of the 
cell size (or 'lattice spacing') — and the inverse volume effects are now sensitive to this new 
scale: The fiducial cell of the homogeneous model is in effect replaced by a physical cell 
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in which the universe can be taken to be homogeneous. However, so far there is neither 
a theoretical principle nor observational guidance on what this scale should be during the 
inflationary epoch when perturbations are generated in these schemes. On the other hand, 
the more striking predictions of these frameworks — such as enhancement of quantum gravity 
signature by several orders of magnitude over the factor H'^/nipi ~ 10~^^ one would expect 
during inflation — appear to depend on the choice of the new scale one makes. Therefore, 
although the underlying idea of lattice refinement is attractive, at present there appears to 
be an inherent ambiguity in the size and importance of the inverse volume effects in this 
setting. Overcoming these limitations is an interesting prospect for future work. 

In these investigations the focus was on quantum gravity effects during inflation. Some 
of the more recent investigations have recognized that, because the energy density is ~ 
10~^^Ppi during this era, quantum gravity corrections during inflation would be too small 
to be observable in the foreseeable future and shifted the emphasis to studying quantum 
gravity corrections from the bounce to the onset of the slow roll. Many of them employ 
a new strategy that goes under the broad theme of 'anomaly cancelation' [47-53]. This 
analysis is based on the Hamiltonian framework and focuses on the constraint algebra. The 
idea is to arrive at the desired LQG theory of cosmological perturbations by studying the 
permissible modifications of the constraints of general relativity that are to encode quantum 
corrections. Recall that two key features of general relativity are: i) the Poisson algebra of 
constraints closes; and, ii) the evolution is generated by these constraints. Keeping these 
features in mind, one proceeds in the following steps: i) One assumes that the 'effective 
theory' that incorporates the LQG corrections would have the same phase space for the 
homogeneous isotropic background and first order perturbations; ii) allows for a modified 
set of constraints on this phase space by making an ansatz for possible modifications; iii) 
calculates the undetermined functions of the background geometry in the ansatz by requiring 
that the constraint algebra should again close; and, finally, iv) defines the desired effective 
dynamics as the Hamiltonian flow of the new, modified constraints. In the quantum theory, if 
the commutator algebra of constraints did not close, there would be an anomaly. Therefore, 
the method is called 'anomaly cancelation' even though one is dealing only with Poisson 
brackets. 

Until recently, the modifications that were arrived at did not change the structure func- 
tions in the constraint algebra of general relativity. In a recent work on scalar perturbations 
[53], on the other hand, the structure function in the Poisson bracket between two Hamilto- 
nian constraint is modified by a function of the ratio p/2pmax where p is the matter density 
in the background and pmax the maximum density in FLRW LQG. This has been interpreted 
to mean that there is a change in the space-time signature at the end of the superinflation 
phase of the background, so that the signature is Euclidean to the past of this event [54]. 

The central idea in this 'anomaly cancelation' strategy is potentially powerful in con- 
straining the type of quantum corrections cosmological perturbation theory can inherit from 
any consistent quantum gravity theory. However, its implementation has some puzzling 
aspects. 

First, the form of permissible modifications of the constraint functionals is chosen pri- 
marily for mathematical convenience and not derived systematically from general physical 
principles. Second, the phase space of the quantum corrected theory is assumed to be the 
same as that in classical general relativity while typically, quantum corrections add higher 
derivative terms to the action which significantly enlarge the phase space. Third, since the 
phase space is kept unchanged, the general analysis due to Hojman, Kuchar and Teitelboim 
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[55] is applicable and their results bring out a puzzling feature of this framework. Hojman 
et al began with 4-dimensional, globally hyperbolic space-times (M, Qab) with non-spinorial 
matter and, using embeddings of a 3-manifold M as a Cauchy surface in M, constructed a 
'hypersurface deformation' algebra V on the space S of embeddings which encodes space- 
time covariance. This is a purely geometric construction without any reference to field 
equations. Then they showed that, on the standard phase space of general relativity (cou- 
pled to the matter under consideration) based on M, there is only one way to represent 
the algebra T> by canonical transformations in a time reversible manner: the representa- 
tion given by the Poisson algebra of the standard Einstein constraint functionals. For the 
anomaly cancelation program, this implies that if the modification of constraint functionals 
is genuine (i.e., not just a field re-definition), then the modified Hamiltonian theory will not 
have a consistent space-time interpretation. Therefore would not be possible to associate 
a well-defined space-time metric to (portions of) dynamical trajectories in this modified 
Hamiltonian theory, let alone examine its signature. Fourth, even if one were to ignore this 
point, signature change is such a drastic effect that it seems difficult to justify the validity 
of the first order cosmological perturbation theory in the subsequent treatment. The fifth 
and a more 'global' limitation arises from the fact that the conceptual underpinning of the 
overall strategy is rather unclear in some of the recent applications. Effective equations are 
meant to incorporate all quantum corrections. Therefore, one would have thought that the 
dynamical equations derived from them would already contain quantum effects. Yet, in some 
works, these fields are quantized again to obtain the power spectrum for scalar, vector and 
tensor modes, as in the standard treatment of cosmological perturbations on classical FLRW 
backgrounds in general relativity. The overall logic underlying this scheme is thus rather 
puzzling. Indeed, already in the homogeneous sector the logical procedure is the opposite: 
one first obtains the quantum evolution equations and then derives effective equations from 
them using dynamics of appropriate, sharply peaked coherent states. 

To summarize, investigations to date have provided useful mathematical infrastructure 
(see, e.g., [47]) and in some cases also the much needed qualitative insights into mechanisms 
through which quantum gravity effects could provide corrections to the standard inflationary 
scenario (see, e.g., [51-53, 56]). The viewpoint that appears to have emerged from the 
'anomaly cancelation' strategy — namely, quantum effects could be neatly encoded in the 
Hamiltonian framework but would lead to a fuzziness if the phase space trajectories are 
interpreted as 4-dimensional space-time geometries — could well be an imprint of some deep 
result on the nature of space-time in LQG. However, it seems fair to say that, at the current 
level of understanding, this and other strategies used so far also have a number of puzzling 
features and face conceptual limitations in their treatment of inhomogeneous perturbations. 

B. Strategy used in this paper 

In this sub-section we will outline the avenue used in this paper to extend the standard 
theory of cosmological perturbations all the way to the Planck regime. As explained in 
section I, the main idea is simply to use the strategy that has driven LQG so far: Construct 
the Hamiltonian framework of the sector of general relativity of interest and then pass to 
the quantum theory using quantum geometry that underlies LQG. In order to bring out 
differences from other approaches, and to address questions regarding truncation and gauge 
choices that are sometimes raised, we will now spell out this strategy using illustrations 
where is has already been successfully used. 
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The FLRW models provide the simplest illustration. Here, one starts out with the full 
phase space F of general relativity (coupled to matter) in the connection dynamics framework 
[57, 58] and truncates it to its homogeneous, isotropic sector Fhi- By fixing gauge, one 
coordinatizes the gravitational part of Fhi simply with a pair (z/, b) of real numbers, where (in 
appropriate units) v denotes the physical volume of the universe (or, if the spatial topology 
is non-compact, of a fiducial cell) and b is its conjugate momentum. This gauge fixing leads 
us to the classically reduced phase space with respect to the Gauss and the diffeomorphism 
constraints. Therefore the reduced phase space carries only the Hamiltonian constraint, 
expressed in terms of just z/, b and matter variables. This provides the starting point for 
quantization a la LQG. In the first step, a specific kinetic framework [59] can be selected 
by a uniqueness theorem [60] along the lines of those in full LQG [61, 62]. By representing 
curvature in the Hamiltonian constraint by holonomies around plaquettes selected by the 
underlying quantum geometry, one constructs the quantum Hamiltonian constraint operator 
[13, 31]. In this construction there are, as is usual in any quantization procedure, some factor 
ordering ambiguities. (Compare, e.g., the Hamiltonian constraint in [13] (which is geared 
to proper time) with that in [14] (geared to harmonic time). See also [63].) However, 
these affect only details of the quantum evolution; general features are robust. Finally, the 
effective equations that encode leading quantum effects are systematically derived starting 
from the quantum theory [40]. 

In the generalization to the homogeneous but anisotropic (i.e., Bianchi) models, one again 
follows the same conceptual procedure: appropriate truncation of the phase space F to Fhom, 
eliminating the Gauss and diffeomorphism constraints by passing to the classically reduced 
phase space, and passing to quantum theory via LQG techniques [21, 23, 24]. The singularity 
resolution persists. Furthermore, there is a detailed consistency check: by tracing over the 
anisotropy degrees of freedom of the Bianchi I model one recovers the quantum Hamiltonian 
constraint of the FLRW model [21]. As one would expect, in presence of anisotropies, 
the quantum dynamics is significantly richer [64-66]. Note that the starting point in this 
analysis is again a truncated phase space Fhom of general relativity, but the truncation is 
now enlarged to incorporate physics of interest, namely anisotropies. 

The last example, Gowdy models, [25-29] illustrate our strategy most closely because 
now the truncated phase space Fcowdy is an infinite dimensional subspace of F: the model 
allows for inhomogeneities induced by a class of (fully non-linear) gravitational waves, in 
addition to those in the matter fields. The constraints are now rewritten in terms of an 
infinite number of conveniently chosen 'modes' of canonical variables coordinatizing Fcowdy, 
obtained by appropriate gauge fixing. Thanks to this choice, we are led to a classical 
reduction of the phase space with respect to the Gauss and the purely inhomogeneous parts 
of the diffeomorphism and Hamiltonian constraints. Thus, one is left with only 'global' 
Hamiltonian and diffeomorphism constraints (corresponding to a homogeneous lapse and 
shift). One then passes to quantum theory using a 'hybrid' scheme where one employs 
LQG quantum kinematics for homogeneous modes and a Fock-type quantum kinematics for 
the inhomogeneous modes representing gravitational waves. This is an internally consistent 
quantization. One finds that, thanks to the quantum geometry effects, the singularity of 
general relativity is resolved. This is in striking contrast to the early attempts predating 
LQG, where singularity could not be resolved (see, e.g. [67-70]). It is interesting to note 
that, if one 'switches off inhomogeneities, Gowdy models reduce to Bianchi I models. Using 
effective equations, it has been shown that the general dynamical behavior of the Bianchi I 
models — including the bounces — carries over to the Gowdy models [26, 29]. This analysis 
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has also provided valuable information on the changes in the amplitudes of gravitational 
waves resulting from the bounce. Finally note that, although one obtains the Bianchi I 
model by switching off inhomogeneities, in contrast to some of the strategies summarized in 
section II A, one does not begin with the inhomogeneous modes propagating on an effective, 
quantum corrected Bianchi I background and then quantize them. Rather, one truncates 
the full phase space F to Fcowdy and quantizes the full Gowdy model, which includes both 
homogeneous and inhomogeneous modes. 

To develop a quantum gravity theory of cosmological perturbations, then, we will continue 
along the same path that has been successfully used so far. Our first task is to identify 
the appropriate truncation of the classical phase space F. By introducing a fiducial flat 
triad for mathematical convenience, we can decompose fields into Fourier modes. With 
this convenient coordinatization, constraints of the full theory can be expressed in terms of 
modes. We will start with a homogeneous background and expand out the deviations from 
it as a sum of first, second, ... nth, ... order terms in inhomogeneities. The truncation will 
now consist of keeping terms which depend on the background and are at most quadratic 
in the first order perturbations. In general relativity, this truncation enables one to study 
dynamics of the background homogeneous space-time and that of linearized perturbations 
propagating on these backgrounds. Back-reaction of these first order perturbations on the 
background (which is encoded in the second order perturbations) is neglected.^ 

The idea again is to use LQG techniques to pass to the quantum theory, using a hybrid 
scheme that broadly mimics the one used in the Gowdy models. (As we will see in subsequent 
sections, some differences arise because, unlike in the Gowdy model, we are now dealing 
with linear perturbations.) It will again be possible to interpret the homogeneous quantum 
Hamiltonian constraint as providing 'evolution' of the quantum state of the background 
geometry in the relational time variable — the scalar field. Solutions to this equation provide 
background quantum geometries on which quantum perturbations evolve. In this conceptual 
setting, one does not start out with classical perturbations on an effective, smooth FLRW 
space-time and then quantize them. The passage to quantum theory is carried out in one 
go for the full truncated phase space that includes both the background and perturbations. 
Finally, as emphasized in section I, much of the technical discussion will be devoted to finding 
a criterion to test whether the final theory admits self-consistent solutions that justify the 
viability of the underlying truncation scheme. 

This overall strategy was briefly reported in [71, 72]. At the same time the 'hybrid 
approach' was used in [73] to study cosmological perturbations in the k=l FLRW context. 
The main conceptual difference between that approach and ours is that we are able to go 
beyond formal considerations in the quantum theory by exploiting the relation between 
quantum fields on quantum geometry and those on a dressed effective geometry discussed 
in section IV B. 

III. TRUNCATED HAMILTONIAN FRAMEWORK 

In cosmology of the very early universe, one generally restricts oneself to the sector of 
general relativity consisting of homogeneous FLRW space-times together with linear pertur- 



^ There are some subtleties in this procedure that arc sometimes overlooked. For a discussion in a simpler 
example, see Appendix A. 
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bations thereon. In much of the cosmology hterature one works with the solution space of 
this truncated theory. However, as pointed out in Ref. [74], the task of finding gauge invari- 
ant variables is more stream-lined in the Hamiltonian framework. More importantly, since 
we are now interested in treating the background geometry quantum mechanically, a natural 
avenue is to follow the Dirac quantization procedure based on phase space. Therefore, in this 
section we will first construct the truncated phase space and then discuss dynamics thereon. 
This will provide a natural starting point to apply the well established LQG techniques in 
the next section. 



Let us begin with general relativity coupled to a scalar field on a space-time manifold 
M = M X M, where M is topologically T"^. For completeness and continuity with the LQG lit- 
erature, we will begin with the connection variables [57] and then pass to the Arnowitt Deser 
Misner (ADM) variables for perturbations that are more commonly used in the cosmological 
literature. 

Let us first focus on geometry. Let qab denote positive definite metrics on M, e°; and 
w*, orthonormal frames and co- frames with respect to Qab, and let Kab denote the extrinsic 
curvature on M. In connection dynamics, the canonically conjugate pair consists of real 
SU(2) connections and su(2) valued vector densities E°- of weight 1, both defined on 
M. (Here indices a,b,c, . . . refer to the tangent space of M and i,j, k, ... to the Lie algebra 
su(2)). They are related to the ADM variables via: 



where q denotes the determinant of Qab, L^, the spin connection determined by e°, A"* = 
KabC^ and 7, the Barbero-Immirzi parameter of LQG. Since we are interested in the spatially 
fiat FLRW background geometries, as usual it is convenient to introduce some fiducial 
structures. Fix a fiat metric Qab, an orthonormal frame e", and the corresponding co-frame 

on M. We will denote by q the determinant of qab and assume that each of the circles in 
M has length i with respect to qab so that the fiducial volume of M is = The natural 
'Gartesian' coordinates defined on M by jab will be denote by x = {xi,X2,X3). 

Even though we are still considering full general relativity, it is convenient to decompose 
the basic canonically conjugate fields into purely homogeneous and purely inhomogeneous 
parts: 



(Here and in what follows, unless otherwise specified, the integrals are over M and dv denotes 
the volume element on it with respect to jab-) Thus, the geometrical part of the phase space 
is naturally coordinatized by quadruples {c,p; a^(x), e°(x)) where c,p are real numbers and 
O-a, e° are purely inhomogeneous fields. The matter field can also be decomposed into purely 
homogeneous and purely inhomogeneous parts: 



A. The Phase space 



Et = ^et and = + 7/^^ 



(3.1) 




/ aie] dv = 0, and / e^w^ dv = . (3.3) 




and 




(3.4) 
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The symplectic structure on the total phase space T is given by 

+ Si(f)62P(^) - 5205iP(^) + / {Siip62TT - S2Tt6i7i) dv (3.5) 

where 6 = {6c, 6p, 5a^, 5e"; 6(j), 5p(^^), 6ip, Sir) denote tangent vectors to T and k = SnG. Thus, 
the only non-zero Poisson brackets between the basic variables are: 

{c, P} = y, {aL(fi), c%X2)} = 6]6'J{x„X2), 
{0, Pw) = 1> Wi^i), tt{x2)} = S{xi, X2), (3.6) 

where 6{xi,X2) = ((l/-\/g)(5(xi, X2) — (1/^'^)) is the Dirac delta distribution on the space of 
purely inhomogeneous fields.^ Thus, the total phase space (F, fl) has a product structure: 

r rhom ^ Tinh) ^ ^hom ~l~ ^inh- 

There are three sets of first class constraints: 

G[A] = jA%x)Gi{x)d^x = jA'{x)DaE^{x)d^x 

Y[N] = jN''{x)Va{x)d^x, S[N] = jN{x)S{x)d^x. (3.7) 

Here, the vector and the scalar constraints are smeared by a shift N and a lapse N and 
the Gauss constraint — which we have explicitly written out because it does not feature in 
the ADM framework — by a generator A* of su(2). The Gauss and the vector constraints 
generate, respectively, internal SU(2) rotations and spatial diffeomorphism, both of which 
are generally regarded as 'kinematic motions.' Dynamics is generated by the Hamiltonian 
constraint. 

Remark: If the spatial topology is rather than T^, the symplectic structure on the 
homogeneous subspace induced by that of full general relativity diverges. Therefore, to 
obtain a consistent phase space description, one has to introduce a cubical fiducial cell C 
aligned with these axes and with edge-length i and restrict integrations of homogeneous 
fields to it. This is an infrared cut-off, to be removed, as usual, in the final results (see e.g. 
[31]). 



B. Expansions around the FLRW subspace 

In physical cosmology one is primarily interested in a neighborhood of the 4-dimensional, 
homogeneous sub-space Fhom of the full phase space F. To fix notation for the rest of the 
paper, we begin with certain expansions of fields near homogeneity. While this technique is 
well known in the perturbation theory around black holes, it appears not to be as familiar 
in the cosmology literature. 



^ The 1 factor ensures that the Poisson brackets are compatible with the fact that the perturbations are 
purely inhomogeneous, i.e., satisfy (3.3). 
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Consider curves 7[e] in F parameterized by e, with e g] — 1, 1[, say, which pass through 
Fhom at e = 0: 



A\[e]{x) = d-'ujl + ea^'^\{x) + ■■■ + -, a^'^'ai^) + • ■ ■ 

n 

EtMix) = VI [Pr'et + e e(^)^(f) + . . . + ;^e("),"(f) + . . .] 

n 

$[e](x) = + e(^(i)(x) + ... + -(^(")(f) + ... 

n[e](f) = [r3p(^) + e7r(i)(f) + . . . + -7r(")(f) + . . .] (3.8) 



where a^^^l^, c^^^i, ff^^K tt^^^ are purely inhomogeneous tangent vectors to the curves 7[e] point- 
ing away from Fhom- Here e is a mathematical parameter that keeps track of the 'order' of 
the perturbation (it could be taken to be y/R, but for simplicity we will assume it to be 
dimensionless here). Geometrically, keeping only the first order perturbations corresponds 
to considering the normal bundle over the homogeneous subspace Fhom of T (since purely 
inhomogeneous tangent vectors are 'orthogonal to' Fhom in the L^-norm). Retaining terms 
only up to the nth order corresponds to considering the nth order, inhomogeneous jet bun- 
dle on Fhom- (Appendix A summarizes the meaning and utility of this expansion procedure 
using the simpler example of the theory.) 

Of special interest are the curves 7[e] that lie in the constraint hypersurface of F. We will 
use a collective label C{A, E,^,]l) for the smeared constraint functions (suppressing the 
smearing fields for simplicity). Then along each curve 7[e] the constraints become C[e] = 0. 
By Taylor expanding in e we obtain a hierarchy of equations: 

CU=o = 0, ^l.=o = 0, ... ^le=o = 0, ... (3.9) 

The zeroth-order equation in the hierarchy, C\e=o = is just the restriction of the full 
constraint to the homogeneous subspace Fhom- As noted above, because of gauge fixing, the 
Gauss and the vector constraints vanish identically and we are left with only one non-trivial 
constraint smeared with a homogeneous lapse [31]: 

2 

SolAThom] := N^om[-^c'\p\'^ + ^] = 0. (3.10) 

i^T 2\p\2 

(As noted in section I, in this paper we have set the potential V{^) of the scalar field 
to zero. See the Appendix in [2] for inclusion of the potential). As usual, although the 
lapse is homogeneous, it can depend on dynamical variables; see section HID for a further 
discussion. Eq (3.10) is a non-linear but algebraic equation constraining the homogeneous 
fields. 

The first order equation, [dC / de]\e=o = 0, is a linear partial differential equation (PDE) 
for a^^^a{x),t^^^i{x), Lp^^\x) , 7i^^\x) . For example, the Gauss constraint yields: 

/Af„h [id'^^^'Kb + cer" e(i)„„ -pr' h'^'^a^'^n] dv = (3.11) 

where we have converted the internal indices i,j, . . . into tangent space indices a,b . . . using 
the fiducial frame and co-frame and indices are raised and lowered with the fiducial 
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metric qab- Only purely inhomogeneous smearing fields contribute because the first order 
perturbations are all purely inhomogeneous. This equation involves aS^\b-i ^^^\b linearly but 
also contains the background variables c, p which solve the zeroth order constraint (3.10). 
In this equation (c, p) happen to enter linearly. But in general the coefficients can be 
complicated non-linear functions of the background fields. Thus, for example, quadratic 
combinations appear in the first order vector constraint, 

/ ^^nh [ip{S^a^^\b - 9aa(i)) - cper"a(i)„„ - c^fe^^^e^^)™,, - K7P(<^)9,(^(i)] dv = (3.12) 

where a = dab^"'^ is the trace of aab- Note that the first order constraints are linear 'homo- 
geneous' PDEs in the sense that there are no (zeroth order) source terms on the right hand 
side. The structure of the first order scalar constraint is the same. 

In the second order equations on the other hand are 'inhomogeneous' as PDEs since the 
zeroth and first order fields now act as 'sources'. These equations determine the 'Coulombic' 
parts of the second order fields in terms of the zeroth and first order ones. For example, the 
second order Gauss constraint is given by 

/A^ [ld-t^''\, + cer" e(2)^„ _ prier"a(2)„„] dv = -f A'^ [iL'^a'^^^h t^'^^] dv . (3.13) 

While this is again a linear PDE — in fact the operator on the left side is the same as that 
in (3.11), but now acts on the second order fields, (a^^^^, e^^^^) — there is now a source term 
which is quadratic in the first order fields which are already known as solutions to (3.11)."^ 
The structure of the second order vector and the scalar constraints is the same. Finally, 
note that the nth order equations in this hierarchy constrain only the nth order fields; they 
do not impose further conditions on lower order fields. This pattern continues to all orders. 
In particular, to obtain the full set of constraints on the zeroth and the first order fields, we 
need to solve only C\e=o = and [dC/d(f]\e=o = 0. This property will be important in what 
follows. 

The key point about the hierarchy is that it greatly simplifies the problem of solving the 
complicated, non-linear PDEs C{A, E;^,]!) = 0. For n = we obtain a non-linear but 
algebraic equation. For n > 0, each equation in the hierarchy is a linear PDE for the nth 
order fields, with the same linear differential operator on the left side but order dependent 
source terms which are already determined by solutions to the lower order equations in the 
hierarchy. The value of the system lies in the hope that by truncating it to a low order, one 
would obtain a good approximate solution to the full system with small inhomogeneities. 
(For further discussion, see Appendix A.) Although it is not easy to rigorously control the 
approximation, this truncation scheme has proved to be a valuable and indispensable tool 
in cosmological and black hole sectors of classical general relativity. In practice, the domain 
of validity of the chosen truncation is checked by self-consistency: one only verifies that, if 
the truncation is of order n, then the source terms in equations governing (n + l)st order 
fields are negligibly small compared to the fields that are kept. 

What is the relation between this hierarchy of constraints and dynamics? On the full 
phase space (F, Q), dynamics is generated by constraints. In particular, in the homogeneous 



^ Because the source terms on the right side of (3.13) are quadratic in first order perturbations, it follows that 
the integral on the right side does not necessarily vanish if the smearing field A" is purely homogeneous. 
The left side of (3.13) now implies that the second (and higher order terms) in our e-expansion (3.8) 
cannot be assumed to be purely inhomogeneous. 
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sector, the scalar constraint smeared with a homogeneous lapse can be thought of as gener- 
ating 'pure time evolution'. What happens if we truncate the theory at first order? Then, 
the truncated phase space FTrun will be the normal bundle over Fhom- The Hamiltonian 
flow generated by S[A''hom] on T is tangential to its homogeneous subspace Fhom- It suffices 
to consider this flow in an arbitrarily small neighborhood of Fhom- Under this flow, tangent 
vectors v = (a^^^^, e*^""^^", tt^^^) at any point on Fhom also have an unambiguous evolution. 
Thus, given a specific Vo at a point, say 7(^0), of any dynamical trajectory 7(t) on Fhom, 
the trajectory 7(t) can be unambiguously lifted to a trajectory in Fxrun, passing through 
{'j{to),Vo)- We will see in section HID that this lift has a simple geometrical interpretation 
in the phase space. 

Remark: In the LQC literature, in place of the e expansion (3.8) one often uses a de- 
composition of the type A'^ = ci'^u^ + 5^4^, etc. Then the 'perturbations' SA"^, . . . include 
terms of all orders n > 1 in our e expansion. Therefore, in light of footnote 4, it is not con- 
sistent to assume that these perturbations SA^^, . . . are purely inhomogeneous. On the other 
hand, if one allows them to have homogeneous parts, then the Poisson brackets between the 
unperturbed and perturbed fields do not all vanish and so the symplectic structure is more 
complicated. This complication is often overlooked. In constraint equations, keeping terms 
linear in these fields is interpreted as the first order truncation, keeping terms quadratic as 
second order truncation, etc. In this scheme, the second and higher order truncation would 
lead to non-linear PDEs on perturbations SA]^, . . .; the simplification that was achieved in 
the e-expansion (3.8) would be lost. Truncations in the two schemes are equivalent only 
to the linear order. These considerations apply also to the Hamiltonian treatment of cos- 
mological perturbations in terms of metric variables, used outside the LQC literature as 
well. 



C. First order truncation 

Since the temperature fluctuations in the cosmic microwave background are only one 
part in 10^, much of the literature on the early universe has focused on FLRW backgrounds 
with first order linear perturbations thereon. We will now analyze this truncation in some 
detail. The relevant phase space is the normal bundle FTrun = Fo x Fi where Fq = Fhom 
is the 4-dimensional homogeneous subspace of the full phase space F and Fi is spanned by 
the first order fields in the expansion (3.8). It turns out that the description of the FLRW 
quantum geometries is easier in terms of variables (b, u), rather than the original (c, p): 

c 4|pP/^ 2 

The geometrical meaning of these variables is as follows^: The physical volume of the universe 
is given by a^i^ = 2TT'~f \u\ ip^, where a is the scale factor, and on any solution b equals 
the Hubble parameter modulo a multiplicative constant, b = 7(a/a). Thus, we will now 
coordinatize Fq by (z/, b; </),p(0)). For simplicity of notation, from now on we will drop the 



^ Like p, the variable ly takes values in the entire real line, positive values corresponding to triads ef with the 
same orientation as the fiducial e°, and negative values corresponding to triads with opposite orientation. 
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sujfix (1) on the first order perturbations. Thus, 

(l/,b, aab,<^ab, V^,V1") G TTmn • (3.15) 

The perturbations Poisson commute with the background fields and the Poisson brackets 
among themselves are given by (3.6) so that f^xrun = + ^i- Thus, mathematically, 
(rxrun, i^Trun) = (r,fi). But the physical interpretation 0/ a^, e^, vr is different: in FTrun 
they represent only the first order perturbations, i.e, just the coefficients of e in the expansion 
(3.8). As a result, the constraints and dynamical equations they satisfy are very different 
from those of the full theory. 

As noted in section II B, to zeroth order we only have the scalar constraint, smeared with 
a homogeneous lapse. In terms of v, b variables it becomes: 

So[iVhom] = iVhom [-^ bV + ^] = 0. (3.16) 

To first order in e, we obtain linear equations on first order fields with functions of b, u, p(^), 
as coefficients. Let us first focus on the Gauss constraint (3.11): 

/ AU [id'^t^'Kb + ch^^ c(i)^„ -pi-' er"a(i)^„]d{; = . (3.17) 

It generates the following infinitesimal gauge transformations: 

aab aab - daAb + Ci-'eabm^"", ^ab tab + P T'^Lbra^'^ ■ (3.18) 

Thus, the symmetric part of tab is gauge invariant and furthermore it has a simple interpre- 
tation. Since E^^ = qq"-'^ in the full theory, where q°''^ is the physical 3-metric on M and 
q its determinant, it follows that 

H<^b) = (l«6 - 1 where qab = q^b + eqafe + ^(e^) . (3.19) 

(Note that q^^jj = (a^) qab is the zeroth order physical metric. It is purely homogeneous 
while is purely inhomogeneous. q := qab^"^.) Next, recall from (3.1) that A^(e) = 
rjj(e) — 7-R'*(e). Linearizing this equation one finds that the first order Gauss constraint 
serves only to determine the skew symmetric part of 

^ab:=u;bi^K{e)U=o (3.20) 
ae 

in terms of tab and the zeroth order fields. Thus imposition of this constraint implies that 
R(^ab) is the free part of M.ab- Furthermore it is gauge invariant. Hence the reduced phase 
space with respect to the first order Gauss constraint (3.11) is coordinatized simply by the 
linearized metric qab and the linearized extrinsic curvature ^ab, both of which are symmetric 
tensors. The three Gauss constraint have removed three degrees of freedom from each of 
the two first order fields (o^, e"), taking us from the linearized connection variables to the 
linearized ADM ones, (qab^^ab)- This is precisely the structure one expects from the Gauss 
reduction of the full theory [57]. 

It now remains to impose the first order vector and scalar constraints on the linearized 
pairs {qab, ^ab)- For this we can draw on the huge body of existing literature. Let us begin 
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by fixing the notation. The ADM variables are {qah^P""^ = y/(j{Kab — Kqab))- To first order, 
they can be expanded as 

9ab(e) = gi? + eqa6 + . . . and p'^^(e) = (P^"^'^' + ep"' + • • •) (3.21) 

where the zeroth order fields are given by q^^^ = a^qab and p^°^°-^ = —[ab/n'y) q°'^. As in the 
cosmology literature, we can expand out the first order fields using Fourier transforms: 



fce£ fce£ 



Pabix) = J ^Pab{k) e^^- vr(f) = j^l^ n^''-"" (3.22) 

where £ is the lattice defined by G ((27r/^)Z)^, A; 7^ 0, Z being the set of integers. 
(The zero k is excluded because, by construction, the fields are all purely inhomogeneous.) 
Since all four fields in the x space are real, their Fourier transforms satisfy the relations 
^ab{k) = c\*fj{—k), etc. It is convenient to expand out the metric perturbations into scalar, 
vector and tensor modes [75]: 



~qab{k) = Sf qab + Sf {kah - Iqab) + V2Vp k^aXb) 



3 

V2VP kaVb) + ^ Tf^ {XaXb - VaVb) + V^TI'^ {X(ayb)) (3.23) 



and. 



Pab{k) = ^ ) qab + ^ Pf\3kakb ' qab) + V^pf ^ k^^b) 

+ V^pf ) kt^aVb) + ^ pf'\xaXb - VaVb) + V2pf ^ (XiaVb)) ■ (3.24) 

(Here k is a unit vector in the k direction and k, x, y constitutes a field of orthonor- 
mal triads with respect to qab in the momentum space. Throughout, indices are low- 
ered, raised and contracted using qab-) Then, the canonically conjugate pairs are 
(4^\p^y, {Sf,p%), ... ((P,.,vr_g), with Poisson brackets {S^^K P%} = 
These fields are subject to three vector constraints and a scalar constraint, each of which 
is linear in these fields but also contains background fields. As is well known, one can pass 
to the reduced, truncated phase space f by solving them and finding gauge invariant 
variables [74, 76]: 

f5iL = {(%.rf,Tf;pf,pf',pf')}^ (3,25) 

r\f r\f Fxi F\i f\i 



f and Tf^ 

k k 

Q^, the gauge invariant Mukhanov-Sasaki variable, is given by 



Here Ti and are the two tensor modes (which are automatically gauge invariant) and 
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Its conjugate momentum is given by 

_ ^ , ^ - ^^'^'^) _ I^ZML _ ^^'^(^) 1 s^^> (3 27) 

The initial configuration variable qab had six degrees of freedom and the scalar field ip, one. 
Each of the four linearized constraints reduces the configuration degrees of freedom by one, 
leaving us with 3 degrees of freedom {Q,T^^\T^'^^) in the reduced configuration space. 

Let us summarize. Note first that in the passage to the reduced phase space FTrun, we 
used only the kinematical structure on FTrun; dynamical equations were not used. Thus, 
the procedure is completely analogous to that followed in other LQC models to extract the 
physical degrees of freedom (see II B). FTrun has the form 

FTrun = F,xf«, with (j., b; 0, p(^) ) G F„ (Q,^,Tf ,Tf ;p?,pf\pJ^)Gf« (3.28) 
These variables are subject to only the zeroth order scalar constraint: 

So[iVhom] = iVhom [-^ bV + ^] = . (3.29) 

We have already taken care of the first order constraints. As noted above, the second and 
higher order constraints do not restrict the first order variables. 

Remarks: 

1. In the CMB, we can only observe modes up to a (finite) maximum wave length Aq which 
equals the radius of the observable universe at the surface of last scattering. Therefore, it 
is physically appropriate to absorb modes with wave-lengths A > lOAo in the homogeneous 
background. This amounts to putting a physical infrared cut-off on perturbative modes, 
making the arbitrariness in the choice of the radius of the 3-torus irrelevant. In phe- 
nomenological applications, these considerations have to be folded into the calculation of 
the renormalized energy density and checking self-consistency. 

2. We use the Mukhanov-Sasaki variable — rather than the curvature perturbation 
TZj: — to coordinatize the reduced phase space because it is better suited for our discussion of 
inflation in [2]. There, we find that TZj^ is not well-defined along effective trajectories because 
it carries p(^^) in the denominator and vanishes at 'turning points' of the infiaton. on 
the other hand is well defined all along these trajectories. In absence of a potential V{(j)), 
as in this paper, they are related just by a constant: TZj: = ^/k/Q Q^. 



D. Truncated Dynamics 

As mentioned in section III A, to obtain dynamical trajectories on FTrun one has to proceed 
in two steps: First obtain a dynamical trajectory on Fq and then lift it to FTrun- 

On F„, 

dynamics is generated by the scalar constraint S[A'^hom]; i-e., the Hamiltonian vector field is 
just the restriction to Fo of the full Hamiltonian vector field^ X°' = f2°''9/3S[A^hom] on F: 



^ Here the Greek indices refer to the (infinite dimensional) tangent space to the full or truncated phase 
space. They are to be regarded as abstract indices a la Penrose. This index notation can be avoided; it 
is used only as a pedagogical aid. 
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X"|r„ = flo'^ di3E>o[Niiom]- Full dynamics on T unambiguously induces a flow on fxrun- The 
general procedure for non-linear systems, summarized in Appendix A, now implies that the 
time evolution on the truncated phase space FTrun is given by the dynamical vector fleld 

Xgy, = nf dpSolN^om] + fif a^§'2[^hom] , (3.30) 

where §2 is obtained from the coefficient S2[A'^hom] of in the expansion of the full scalar 
constraint S[A''hom] on F, by keeping just those terms which are quadratic in the first order 
quantities (ai^''*, e'^^^'^j, f/?*^^-', vr^) and discarding terms linear in the second order quantities 

It is important to note that it is only E>2 [A'^hom] that is constrained to 
vanish; there is no constraint on S2[A'hom]- Finally, because §2 depends on the background 
fields, A'pyj^ is not generated by §0 + ^2; i.e., -^Dyn 813(810 + ^2) ■ But one might wonder 

if one can find another 'effective Hamiltonian' that generates this dynamical flow. The 
answer is in the negative: One can check that X^^^^ does not Lie drag the total symplectic 
structure Q = Qo + ^1 on the truncated phase space, whence it is impossible to flnd an 
effective Hamiltonian. 

Remark: We have spelled out these results because there has been some confusion 
in the recent LQC literature. Quantum dynamics is often assumed to be captured by a 
quantum constraint Cfj"^ = 0. Our detailed discussion of truncation and dynamics of the 
truncated system implies that this assumption cannot be justifled. First, as Eq. (3.30) 
shows, the dynamical vector fleld ^Dyn Fxrun is not generated by a constraint (or indeed 
by any Hamiltonian). Therefore there is no reason to expect that the correct quantum 
dynamics could be recovered by imposing any quantum constraint. Second, even the part 
Q"^^ di3E>2[N^om] of X^y^ describing the evolution of perturbations is not generated by the 
second order constraint S2[A'hom] but by S2[^hom], which is unconstrained. 

Finally, we are interested only in the dynamics of gauge invariant variables, i.e., only in 
dynamics on Fttuii. As one would expect, the_functions So and S'2 that determine X" are 
gauge invariant and therefore projects down to Fxrun. In particular, S2 = §2^^^+^2^"^^''+^2''"^^'' 
and the form of the last two terms is exactly the same. Therefore, from now on, we will 
denote the two tensor modes collectively by T. Then [74], 

k 
k 

Note that if we were to replace l^/HQk by Th (and the corresponding momenta by a reciprocal 
factor to maintain the Poisson brackets), §2^^'' [A^hom] reduces to §2^^'' [-^hom]. Therefore, it 
suffices to focus just on ( one of the two tensor modes ) 7^ and its conjugate momentum which 
we will denote simply by p^. 

The resulting dynamical evolution can be best understood in geometric terms as follows. 
Note flrst that FTrun is naturally a bundle over the homogeneous phase space To- Fix an 
integral curve of X"^ on the 'base space' Fq, i.e. a homogeneous solution. Fix a point 
7o = (z/°,b°; 4>"yP(^)) on this trajectory and a tangent vector 7^°,p| at that point. To 
describe the evolution of this perturbation on the chosen background trajectory, we need 
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to lift the trajectory from To to Fxrun- This can be done simply by considering the integral 
curve of X^^^ passing through the point (z/°, b°; (t)°-,p°^y 7^°, p|) of FTrun- 

Next, let us list the commonly used lapse functions iVhom and the corresponding time 
variables: 

• ^hom = 1 corresponds to proper time t so that the physical space-time metric has the 
form ds^ = —dt^ + a?dx^. 

• -^hom = ci corresponds to conformal time rj so that the physical space-time metric has 
the form ds^ = o?{—drf -\- dx^). This is the most common choice in the cosmology 
literature. 

• -^hom = 0? corresponds to harmonic time r which satisfies the wave equation, Dr = 
0. The physical space-time metric now assumes the form ds^ = —a^dr^ + o^dx"^. 
The zeroth order scalar constraint So[^hom] takes the simplest form with this choice. 
Therefore, harmonic time is commonly used in LQC not only for the /c = 0, A = FLRW 
cosmology but also for models that admit spatial curvature, a non-zero cosmological 
constant and anisotropies [14, 21, 23, 24, 31]. 

• N = (V^/p(0))a^ corresponds to choosing the scalar field itself as time. It turns 
out that irrespective of the initial choice of lapse in the classical theory, the quantum 
scalar constraint has a form that naturally leads one to use as a relational or internal 
time variable [13, 14]. In the case now under consideration, where the scalar field (j) 
also satisfies the wave equation 00 = 0, the internal time defined by and the 
harmonic time r are related just by a constant in any classical solution: = {P{4,)/V)t. 
However, since the constant varies from one solution to another, in the quantum theory, 
conceptually, the two choices are quite different. In LQC it is simplest to begin with 
the harmonic time in the classical theory and reinterpret the quantum scalar constraint 
as providing time evolution in the relational time 0. 

Since one generally uses the conformal time r] in the cosmology literature, we will conclude 
by writing down the equation of motion that follows from (3.31) with the choice A^hom = Ci'- 

7^" + 2 ^ 7^' + er-^ = (3.32) 

where a 'prime' denotes derivative with respect to the conformal time rj. One often rescales 
Tj: to obtain a field Xk "with physical dimensions of a scalar field, 

X,.:=^r,, (3.33) 

for which, furthermore, the equation of motion resembles that of a harmonic oscillator (with 
time dependent frequency): 

X'^ + (-T + - 4 + + ' (3.34) 
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where R denotes the scalar curvature of the background homogeneous space-time/ This 
form of the perturbation equations has also been exploited in [77-79] in singling out a 
physically motivated quantization. 

Remark: In our presentation we began with the connection variables employed in the 
LQG literature because conceptually they are essential for the new quantum kinematics 
used in LQG, and the subsequent treatment of the singularity-free quantum dynamics of 
our background FLRW space-times. But we passed to the ADM variables for first order 
perturbations by solving the Gauss constraint because much of the cosmological perturba- 
tion analysis is carried out in terms of these variables. For a treatment of cosmological 
perturbations in connection variables, see in particular [76, 80-84]. 



We will now use the phase space Fxrun as the point of departure for quantization. Since the 
truncated second order constraints for the scalar and two tensor modes are identical (except 
for unimportant numerical factors), as in section III D we will focus just on one tensor mode 
7^. Thus, from now on, Fivun will be taken to be spanned by three canonically conjugate pairs 
(z/, 6; (p,P(if,); 7^, pg), the first two representing the background, and the third representing 
the first order perturbation. They are subject to the single constraint So[A'hom] = (see Eq 
(3.29)). Dynamics is governed by the vector field ^oyn = ^^^dpEo + r2"'^5/3§2 where §2 is 
given by (3.31). 

In this section we present the general program without entering into details of how the 
states and operators related to first order perturbations are defined. In particular, the 
expression of the Hamiltonian operator dictating the dynamics of first order perturbations 
is formal. The precise definitions of these states, operators and the necessary regularization 
procedure are provided in the next two sections 



Let us begin by recalling the quantum theory of the homogeneous sector (for details, see 
e.g. [13, 14, 31]). As mentioned in section II, LQG provides a kinematical Hilbert space "Hkm 
which, as in full LQG [61, 62] is uniquely determined by the physical requirement indepen- 
dence w.r.t. background (and fiducial) structures [60]. In the configuration representation, 
kinematical quantum states are given by wave functions \l/(i^, 0). As noted in section HID, 
dynamics is simplest if one uses harmonic time by setting N = a^. Then, the scalar con- 
straint becomes a well-defined operator on 7/kin- Physical states "^oi^^, 4>) are annihilated by 
this constraint, i.e., they satisfy 



e^o(z/, 0) = -— z/[(z/ + 2A)*,(i/ + 4A, 0) - 2z/*o(i/, 0) + (z/ - 2A)^,(z/ - 4A, 0)] . (4.2) 



IV. QUANTUM THEORY: MAIN STEPS 



A. Background quantum geometry 




(4.1) 



where the action of 6 is given by 



For the scalar perturbation, -0^: — aQj: has the same properties as Xk'- V'fe 3'lso has the physical dimensions 
of a scalar field and satisfies the same equation of motion as (3.34). 
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Thus, is a second order difference operator that acts only on the argument v of \l/o; with 
step size 4A where = 4^371741 is the 'area gap' of LQG. 6 is self-adjoint and positive 
definite on "Hkin [85]. As is common in the Dirac quantization procedure, none of the solutions 
to (4.1) are normalizable in the kinematic inner product. But there is a standard 'group 
averaging' method to endow the space of solutions with a physical inner product [86-88]. 
The resulting physical states '^o £ ^phy ^^^'^ out to be solutions to 

-ifid^^iy, (p) = Ho^o{i^, <P) where Ho = hVO, 
which are symmetric under z/ — t- — z/ and have finite norm 

ll^o|P = - V |H"M*o(i^,0o)P, (4.4) 

where (po is any fixed instant of the internal time 0. (The scalar product defined by (4.4) 
is insensitive to the choice of 0o-) Heuristically (4.3) can be thought of as the positive 
frequency square- root of (4.1) in the internal time 0. 

The most noteworthy feature of this outcome is that the quantum Hamiltonian constraint 
is naturally de-parameterized: its form suggests that the scalar field can be interpreted 
as a 'relational or internal time' with respect to which physical states '^o evolve. Thus, as 
one would hope, the imposition of the quantum constraint a la Dirac has naturally led us 
to dynamics. Interestingly, while our use of a lapse function corresponding to the harmonic 
time r simplified the form of Ho-, it was not essential to arrive at this interpretation of 0. 
Indeed this interpretation was initially arrived at with lapse set to 1 corresponding to proper 
time [13] (but a more complicated form of Ho-) 

Thus, as is usual, the Dirac quantization procedure has naturally led to a Schrodinger 
picture in which the scalar field is simply a time parameter and the sole dynamical variable 
is V that determines the volume of the universe via 

v^!o{y) = 2'K^^lM'^o{y)- (4.5) 

In the Heisenberg picture, the volume evolves in time from the bounce via 

y(0) = (Ai/mo{4>-M (27r74i|z/|) e-(*/^)^°('^-^B) 

Every element \&o of ^phy represents a 4-dimensional quantum geometry. However, for 
our purposes, only a subset of these states is relevant. Choose a classical, expanding FLRW 
space-time in which ^ h (in the classical units G=c=l) and a homogeneous slice at a 
late time = 0o, when the matter density and curvature are negligibly small compared to 
the Planck scale. This defines a point 70 in To- Then, in the Schrodinger representation, one 
can introduce 'coherent states' \E'o(i^, 0o) in "Hp^y which are sharply peaked at 70 [12]. By a 
quantum background geometry, we will mean a physical state ^ o{i^i 0) obtained by evolving 
these initial states using (4-3). There is a large class of such states and our considerations 
will apply to all of them. One can show that these states remain sharply peaked on the 
classical trajectory passing through 70 for all > 0o. In the backward time-evolution, 
they do so only till the density reaches a few hundredths of the Planck density. Even 
in the deep Planck regime the wave function remains sharply peaked but now the peak 
follows an effective trajectory which undergoes a quantum bounce. At the bounce point the 



(4.3) 



(4.6) 
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matter density attains a maximum, Pmax ~ 0.41ppi. While there is good agreement with 
general relativity once the matter density falls below a few hundredths of the Planck density, 
Einstein's equations break down completely in the Planck regime. But the quantum state 
(and even the effective trajectory) remains well-defined throughout the entire evolution, 
including the Planck scale neighborhood of the bounce [13, 14, 31]. 

Finally, in the Heisenberg picture, we can define the space-time metric operator. Recall 
first that in the classical theory the lapse function corresponding to the scalar field time 
is given by iV = {(j)) p'^^y Since P{(f))'^o = Ho^o for any \&o G "Hp^y, in the Heisenberg 
representation, the metric operator is given by 

gabdx^dx'' = = H~H^ {(j)) R-^ dcj)"^ + 0^(0) dx^ (4.7) 

where we have used a symmetric factor ordering in the first term and defined the (positive 
definite, self-adjoint) scale factor operator in the Heisenberg picture via: 

ra{<P) = [Vm^- (4-8) 

In the Heisenberg picture, the geometry is quantum because the metric coefficients are now 
quantum operators on "Hphy- 



B. Perturbations on the cjUcintiim geometry ^ o 

Because we were able to reinterpret the quantum constraint equation §o^o = as provid- 
ing an evolution of physical states in the internal time variable 0, we were naturally led to 
work in the Schrodinger picture for the homogeneous background geometry. Furthermore, 
because o represents a quantum state of the background geometry and if) of perturbations, 
it is natural to assume that the total state has a simple tensor product structure 

*(i/, Tg, 0o) = *o(i^, 0o) ® ^(7i, 0o) (4.9) 

at some initial time </)o. Then, because the back- reaction is neglected, the evolution of \&o 
is dictated just by (4.3); it is insensitive to the form of if). (This is entirely analogous the 
situation in the classical theory.) Therefore the tensor product structure is preserved under 
evolution. As in the classical theory, our task is to evolve if) on the specified background 
quantum geometry \E'o, i.e., to lift the given homogeneous 'quantum trajectory' \&o(i^, 0) to 
a 'trajectory' \&o(z/, 0) ® ^/;(7g, 0) of the truncated quantum theory, where the background 
state is the given one. To carry out this task we need to first complete two preliminary 
steps. 

On the classical phase space FTrun = Fq x F*^^), the part Vl'^^d^'B^ of the dynamical vector 
field ^Dyn dictates how perturbations propagate on a homogeneous background solution 
v[(f)). Now, §2 depends not only on the perturbations (7^, p^) but also on the time dependent 
scale factor of the background solution. Therefore, to construct the operator Sg? if is simplest 
to work in the 'interaction picture' where the background scale factor operators evolve in 
the relational time (f) and the background state o is frozen at a time, which we will take 
to be the bounce time (p = (p^. The first preliminary step is to carry out this passage to the 
interaction picture via 



*int(i^, Tg, 0) = e-(^/^')^°(<^-^B) (vp„(i., 0) ® ^(r^, 0)) 



(4.10) 
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A second step is needed because now the evolution is with respect to the relational time 0. 
Therefore, we have to choose a specific lapse function, A^hom = a^-^^/P((/>), in the expression 
(3.31) of §2 cind then use an appropriate factor ordering to convert it to an operator. In this 
step, we will use the same factor ordering as in the expression (4.7) of the quantum metric 
operator and make a simplification using the evolution equation P(^)\l/o = — "^^f^^^o = -^o^o- 
These two steps, and the form (3.31) of E>2 lead us to the following evolution equation for 
the total state in the interaction picture: 

= \ X^4«:[i7-ivl;,(,.,0j]® [|p^^|V(r^,0)] 

+ ^[(^o"^aV)^o-^)vI/o(z.,0j]®[|r^|2^(rg,0)]. (4.11) 

Let us take the scalar product of this equation with (which we assume to be normal- 
ized), the result is the required evolution equation for the quantum state ?/'(7^, 0) of the 
perturbation propagating on the quantum background geometry '^o'- 

k 

+ ^(^o-^a^(0)/7o'^)|r^|V(r^,0), (4.12) 

where, by construction, the expectation values of the background geometry operators are 
taken in the given state ^ o of the background quantum geometry. 

Note that, at a fundamental level, i\} now evolves in a probability amplitude \E^o of back- 
ground geometries Qahi rather than a fixed Qah- But (4.12) implies that its evolution is not 
sensitive to all the details of the fluctuations of this quantum geometry; it is sensitive to 

only two 'moments' {H~^) and {Ho ^ a^(0) Ho'^). This is so even though (4.12) is an exact 
consequence of (4.11) with no further approximations. Furthermore, since the back-reaction 
of the perturbation on \E^o is neglected within the test field approximation inherent to our 
truncation scheme, nothing is lost by projecting (4.11) along \l/o in arriving at (4.12). That is, 
within the test field approximation, (4-12) captures the full information about the evolution 
of ip that is contained in the original equation (4-11) ■ 

Next, recall that in the standard cosmology literature one regards the quantum state 
i^iT^, 0) of perturbations as propagating on a classical FLRW metric, specified by the scale 
factor aci(0). In the Schrodinger picture now under consideration, the evolution is given by 

^m,^|^{T^, <P) = \Y^ 4«:(p-p |p,H^^(rg, 0) + ^ [p-^^ftm \f^'i^{T^, 0) (4.13) 

k 

where pi^^) can also be expressed using geometric variables using the constraint equation 
satisfied by the background. Comparing (4.13) with (4.12), we find that the evolution of 
the test perturbation 7^ on the quantum background geometry given by \&o(i^, <P) is indis- 
tinguishable from that of a test perturbation propagating on a smooth FLRW background 

"gabdx^dx^ = d~s^ = -iPi^))-^ f d\(j)) d(j)^ + a(0)2 dx^ (4.14) 
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where 




{Ho') 



1 



-1 



and 



(4.15) 



Again, this is an exact equivalence between our truncated LQG and the theory of quantum 
fields on a smooth FLRW geometry determined by d and p(^^y 

At a technical level, the existence of such a simple relation is quite surprising at first. 
But this result should not be interpreted to mean that the standard quantum theory of 
perturbations on a classical FLRW solution to Einstein's equations holds in the Planck 
regime. It does not! For, the a(0) seen by 7^ is very different from the aci(0) of a classical 
solution; in particular, (5(0), p(0)) do not satisfy Einstein's equations. Indeed, their 
expressions involve H; although cjab is smooth, it incorporates quantum corrections which 
are so large in the Planck regime that they tame the big bang singularity. Furthermore, the 
pair does not even satisfy the effective equations of LQC which track the peak of the wave 
function \l/o(i/, 0) of the background geometry. Certain aspects of quantum fluctuations 
inherent in \&o('^, 4>) are absorbed in these tilde fields. Thus, (jab may be thought of as a 
dressed effective geometry that is relevant for propagation of linear perturbations on the 
full quantum (background) geometry determined by \E'o(^, 0)- In retrospect, from what we 
know in other areas of physics, such a result is not entirely unexpected. For example, light 
propagating in a medium interacts with its atoms but the net effect of these interactions 
can be encoded just in a few parameters such as the refractive index of the medium. In 
our case, the 'medium' is the quantum geometry and the tilde variables d,p(^^) encode the 
interaction between this 'medium' and the perturbations. Only two parameters suffice 
simply because the background quantum geometry is homogeneous and isotropic. Already 
in the anisotropic Bianchi models, (an extension of the discussion of [89] implies that) one 
would need additional parameters characterizing the dressed, effective anisotropies. Finally, 
the encoding is rather sophisticated: prior to the calculation, it would have been impossible 
to guess the precise 'moments' of the fluctuations of geometry that are to capture this 
interaction. 

Remark: Several equations in this sub-section closely resemble those in [38]. However 
the conceptual under-pinning is quite different. The discussion in [38] began by assum- 
ing that quantum dynamics can be obtained by imposing the constraint analogous to 
[So -l- S2] ® '0 = 0. This 'quantum constraint' was then expanded and sub-leading terms 
were discarded to arrive at an equation analogous to (4.11). Consequently it was not 
appreciated that, within the truncation scheme, (4.12) carries the full information about 
the quantum evolution of ijj. As discussed in section HID (see also Appendix A), a more 
careful examination has revealed that the strategy of imposing a quantum constraint cannot 
be justifled. Therefore, we adopted a different route here. We mimicked the strategy from 
the classical theory: Just as -^Dyn provides a lift of the dynamical trajectories on Pq to 
To X Pi, we lifted the 'quantum dynamical trajectories' on 'Hp^y '^o<S)4' on "Hpj^y (gj'H^. 
This route is also more direct in that we did not have to discard any terms to arrive at (4.11). 

Finally, let us translate this result using conformal time 77 commonly used in the literature 
on cosmological perturbations. The dressed, effective metric can be written as: 




(4.16) 
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where 

dv=[fd\<p)]pl^]d<p. (4.17) 

Therefore, in the truncated theory, the exact evolution equation for the quantum perturba- 
tion 7^ on the background quantum geometry is given by 

Tl' + 2^fS + k% = Q (4.18) 

where the prime now denotes a derivative with respect to fj. This mathematical equivalence 
simplifies both conceptual and technical aspects of our analysis considerably because the 
well-developed techniques from quantum field theory in curved space-times can now be 
readily imported into the quantum field theory of perturbations ip on quantum geometries 
\l/o. In section VI, we will use this strategy to define in detail the quantum states for 
perturbations and composite operators that are needed to complete the quantum theory. 

Remarks: 

1. We can make a further simplification through a 'mean field' approximation in which the 
fluctuations are ignored. More precisely, let us first recall [13, 14, 31] that even in the Planck 
regime the state \l/o(i^, 0) is sharply peaked on an effective geometry 

£6 

Qabdz^dx^ = df = d^^d(f)'^ + d^U)dx^ (4.19) 

which agrees with the general relativity solution (for the same value of P(0)) for large a(0) 
but has a in-built bounce at a(0)^ = P^(^)/(2^^Pmax)5 with p^ax = 3/(87rG7^A^) ~ 0.41ppi. In 
terms of our quantum geometry state the scale factor is given just by the expectation 
value a(0) = (a(0)) in this state. Suppose we ignore quantum fluctuations, i.e., use a mean 
field approximation in which the expectation values of powers of a and Hq are replaced by the 
same powers of their expectation values. In this approximation, the quantum perturbation 
7^ would seem to propagate on the effective geometry determined by the pair (a(0),p((^)). 
Now, one knows that there exist background quantum geometries \&o which are very sharply 
peaked on this effective geometry even in the Planck regime. If one uses such a \l/o, the mean 
field approximation is excellent for studying the propagation of perturbations on quantum 
geometries under consideration.^ The exact evolution of 7^, on the other hand, sees the 
more sophisticated, 'dressed' effective geometry determined by (a,p(0)). 

2. If one were to use the mean field approximation, the quantum perturbations would 
satisfy (4.18) with the tilde quantities replaced by the barred quantities that refer to the 
effective LQC solutions. At first sight, it may therefore seem one could have arrived at 
these equations simply by perturbing the effective equations of LQC. This could have been 
a viable interpretation had there been a clear set of effective equations in full LQG to 
perturb around backgrounds satisfying the LQC solutions. But as emphasized in sections 
I and II: i) we do not yet have effective equations in full LQG, and, ii) if one were to 
adopt the naive strategy of considering a set of linearized equations in general relativity 



^ In numerical simulations of the evolution of the quantum state ij: of perturbations, for example, if ^ o is 
chosen appropriately, the numerical errors would be much higher than those introduced by the mean field 
approximation. 
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and simply replacing the background solution to Einstein's equations by a solution to the 
effective equations, one faces a very large ambiguity in the choice of equations with which to 
begin and, furthermore, the set of final equations need not be internally consistent when the 
background does not satisfy Einstein's equations. Our procedure is free of these drawbacks 
because we first constructed the quantum theory and arrived at Eq. (4.18) by showing 
an exact equivalence of fields Q, T propagating on the quantum geometry \E'o and those 
propagating on the geometry determined by cjab- Nowhere did this procedure use effective 
equations of LQC 

V. HILBERT SPACES OF STATES OF % 

k 

In this section we will construct the Hilbert space of quantum states of gauge invariant 
perturbations following two different but complementary avenues. The first is geared to 
mathematical physicists and well adapted to the Hamiltonian framework used in our clas- 
sical considerations. The second follows the route that is more often taken in analyzing 
cosmological perturbations. We show that the two are equivalent. Therefore, to ensure 
conceptual continuity and coherence, one can start with quantization given in section VA 
and then use the framework presented in VB which is better adapted to regularization, 
renormalization and numerical simulations. We also provide the explicit expression of the 
2-point function which makes it manifest that the group of space-translations continues to 
be a symmetry in the quantum theory. 

A. The Weyl algebra and its representations 

The classical phase space Tj- for the tensor modes is spanned by canonically conjugate 
pairs (7^, p^). The corresponding operators 7^, satisfy the canonical commutation rela- 
tions and generate the Heisenberg algebra in the quantum theory. For technical simplicity, 
it is convenient to exponentiate them to obtain the Weyl algebra QU whose representations 
provide the required Hilbert spaces of quantum states. 

One begins with the observation that, with each vector (A,/i) G Ft-, one can associate a 
natural linear combination of smeared configuration and momentum operators: 

F(A,/i) :=fi((r,p), (A,/i)) = i J2f^^% - 

k 

= I d^x {ii{x)f{x) - \{x)p{x)) . (5.1) 
Jm 

The Weyl operators VF(A,/i) are their exponentials: 

iy(A,/i) :=es^(^''^). (5.2) 

It is more convenient to work with these exponentials for two reasons. First, while the 
field operators F are unbounded, the W are unitary and hence bounded operators in any 
representation. Therefore one avoids the awkward issues of specifying operator domains. 
Second, the vector space generated by finite linear combinations of the Weyl operators is 
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automatically closed under the Hermitian-conjugation operation W^{X,^) = W{—\, —fi), 
and, more importantly, under the product: 

WiXi, /ii) WiX2, /i2) = /-i'- yy^^^ ^ ^ _ ^5 3) 

Thus, the vector space has the structure of a T*c-algebra. This is the Weyl algebra Q2T. (It 
can be easily endowed the structure of a C* algebra but this will not be necessary for our 
purposes.) 

To find the representations of 211, it is simplest to use the standard Gel'fand, Naimark, 
Segal (GNS) construction [90]: Given a positive linear function (PLF) on QU, the construction 
provides an explicit Hilbert space "Hi and a representation of elements of QU by concrete 
operators on that "Hi. This representation is cyclic: every state in "Hi arises from the action 
of operators representing elements of 21J on a 'vacuum'. In this representation, the PLF 
turns out to be just the vacuum expectation value functional. A natural strategy for linear 
fields is to first find a complex structure on the phase space that is compatible with the 
symplectic structure thereon and then use the Hermitian inner product provided by the 
resulting Kahler structure to define the required PLF on QU [91, 92]. We will follow this 
conceptual strategy but in a manner that retains close contact with the cosmology literature. 
Therefore, prior knowledge of the complex and Kahler structures will not be necessary to 
follow the construction. 

In linear field theories in Minkowski space, one narrows the selection of the positive linear 
functional by requiring that it (and hence the vacuum state) be Poincare invariant. In the 
present case, it is natural to require that the PLF be invariant under the 3-dimensional 
translational symmetry of the background geometry. Such PLFs can be constructed as 
follows. Choose a set of complex coefficients (cfc, fk) (with k >0) such that 

efc/fc - 4fk = 2« Im (cfc/^t) = t for all k (5.4) 

(For a massless scalar field in Minkowski space, = e"""^*"/ \/2uj and fk = {—iu e"""^*")/ \/2uj 
with u = \k\.) Then, one can extract the 'positive frequency' part of any vector (A,/i) G 
Tj- as follows: 

H := -i {fkh - 4h) (5-5) 

so that 

h = (^ka,: + e^a|, and /ig = /^a^ + f^al . (5.6) 

(Note that while A, /i satisfy the 'reality condition' A| = A_^, /i| = the 'positive 

frequency parts' don't: at ^ a_^.) The required PLF is then simply 

(#(A,/i)) = e-^7 J^l'^sl' (5.7) 

Under the action of a translation x — )■ x + {T on M, we have: T'(x) — T(x + {/), and 
p(x) — )■ p(x + v), whence l^(A(x), yU.(x)) — )■ W{\{x — v),ji{x — v)), or, in the momentum 
space, A^ — )■ e~*'^'^'Ag, /ig — )■ e~^^'^fi^. This trivially implies — )■ e~^^""a^. Hence the PLF 
(5.7) is left invariant and the translation is represented by a unitary transformation on the 
GNS Hilbert space "Hi under which the GNS vacuum is invariant. Thus, each choice of 
coefficients (cfc, fk) satisfying (5.4) leads to a representation of the Weyl algebra in which 
the (GNS) vacuum is invariant under translations. 
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B. Space-time description 

At this stage it is convenient to make the construction more exphcit using the famihar 
expansions in terms of the creation and annihilation operators. This entails going to the 
Heisenberg picture. In terms of the conformal time fj of the dressed effective metric of section 
IV A, the space-time operator is represented as T{x,fj). It satisfies the field equation and 
the canonical commutation relations, and is self-adjoint. These properties can be neatly 
captured by an expansion of the form 

k 

We require that efc(r/) should satisfy (4.18): 

el{v) + 2^e',{v) + ee,{v) = 0, (5.9) 

so that T'(x, fj) satisfies the desired equation of motion. (In relation to the more familiar 
expansion in Minkowski space-time, 6^(77) now plays the role of the positive frequency basis 
functions e~*'^*/v^2cJ.) Next, for each k, the space of solutions to (5.9) is two dimensional 
and one chooses a complex solution 6^(77) satisfying the normalization condition: 

^ {ekiv) etiv) - eliv) eiiv)) = ^ . (5.10) 

This condition needs to be imposed only at some initial instant of time fjo] Eq (5.9) then 
guarantees that it is then automatically satisfied for all fj. With this normalization, if the 
time independent operators and At satisfy the commutation relations 

4'] = ^^'%" (5.11) 

then 

fix, v) and m v) = ^.IY1 (^'fc(^)4^ + ■ (5-12) 

k 

satisfy the required canonical commutation relations at any fixed value of fj. In view of their 
properties, one can interpret as annihilation operators, define the vacuum |0) as the state 
annihilated by all A^, and generate the Fock space "Hi by repeatedly acting on the vacuum 

by creation operators At. It is straightforward to calculate the vacuum expectation value of 
the Weyl operator 

W{\n)\f,^ = J^£/^Fe-^Ipe (5.13) 

at any conformal time fjo. It is given by (5.7) where = 6^(7^0) and fk = [a? e'^.(f/o). 
Thus, our second description of the Hilbert space l-ii adapted to the covariant space-time 
picture, is completely equivalent to the first description, adapted to the Weyl algebra that 
was constructed starting from the phase space. The first description serves to bring out the 
conceptual structure of the quantum theory that emerges from the phase space description 
of section III while the second is better adapted to calculations, e.g. regularization of the 
stress-energy tensor discussed in section VI. 



29 



To summarize, for each choice of solutions 6^(77) satisfying the normahzation condition 
(5.10) we obtain a vacuum state and hence a Fock representation of the canonical com- 
mutation relations (or, of the Weyl algebra). Since these representations are completely 
characterized by their 2-point functions, it is instructive to write them out explicitly. Using 



By inspection the 2-point function is invariant under the action of space-translations. This 
is an independent proof of the translational invariance of the vacuum state, now geared to 
the cosmology literature. 

Remark: While each choice of the family of solutions ek{fj) to (5.9) satisfying the 
normalization condition (5.10) determines a vacuum state (and the associated Hilbert space 
Hi of states), this is a many to one map: The families 6^(77) and e^^'= efc(r/) that differ just 
by a /c-dependent phase factor determine the same vacuum. There is a 1-1 correspondence 
between the equivalence classes {ckif])} of families that differ by such phase factors and 
complex structures J on the phase space Ft- which are compatible with the symplectic 
structure and are invariant under the action of space-translations. Thus there is a 1-1 
correspondence between the vacua |0) we have constructed and complex structures on F7- 
satisfying the two properties listed above. 

A natural question about these representations of the Weyl or the Heisenberg algebra is 
now the following: Does a change of the complex structure J — i.e., the choice of 'generalized 
positive frequency solutions' {efc(f/)} — always result in unitarily equivalent representations? 
As is well-known, in general the answer is in the negative. In the terminology used in the 
cosmology literature, in general the vacuum state selected by any one complex structure may 
contain an infinite number of particles corresponding to another complex structure. A priori 
this would be a key obstacle to extracting physics because different choices would in general 
lead to very different predictions. Furthermore, in a general representation so constructed, 
there is no natural prescription to regulate products of operator-valued distributions, e.g. 
T^{x,fi), and hence to define basic physical operators such as the Hamiltonian of (4.12). 

Fortunately, as we will see in section VI, both these problems can be resolved in one 
stroke by imposing certain regularity requirements on the basis functions {e-k{fj)}. Then, 
the representations of the Weyl algebra that result from equivalence classes {ek{fj)} of any of 
the regular basis functions — or 'regular' complex structures J — will turn out to be unitarily 
equivalent. In this sense there is a unique class of unitarily equivalent representations and one 
can work with a unique Hilbert space "Hi. Therefore the more general framework of algebraic 
quantum field theory is not essential in the cosmological context under consideration.^ 



^ See also [77-79] where one arrives at a rigorous uniqueness result using the form (3.34) of equations of 
motion. However, the 'vacua' they are led to consider are not necessarily of 4th adiabatic order whence 
it would be difficult to regularize and renormalize physically interesting composite operators, such as the 
energy density, in that framework. It would be interesting to see in detail the precise relation between 
that approach and the adiabatic treatment pursued here and in much of the cosmological literature. 




(0 I r(fi,r^i)r(x2,?72) I 0) = Y.^^-^''~'"^ eu{vi)el{fi,) . 



(5.14) 



k 
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All the regular translationally invariant vacua and states containing a finite number of 
excitations over any of them belong to the Hilbert space "Hi. However, the translationally 
invariant vacua span an infinite dimensional subspace of "Hi and none of them is preferred 
from a full space-time perspective. Thus, T-Li will not admit a 'canonical' vacuum we are 
used to in Minkowski or (strictly) stationary space-times. 



VI. REGULARITY CONDITIONS ON STATES AND OPERATORS 

In this section we will first summarize the notion of regularity conditions on states ip of 
quantum perturbations and then use it to regulate products of operator- valued distributions, 
such as the ones that appear in the quantum stress-energy tensor or Hamiltonian. There 
exist several methods of regularization. We will work with the adiabatic scheme because it is 
particularly suitable to perform explicit computations, including numerical implementations, 
and can be directly extended to our quantum field theory on quantum FLRW geometries. 
In this paper of course we will focus on the k = 0, A = case but our considerations will 
extend to other contexts such as the k = l FLRW case where the underlying isometrics make 
a mode decomposition naturally available. 

Much of the discussion in the first two sub-sections is taken from the rich literature 
on adiabatic regularization in cosmology (see, e.g., [8, 9, 95, 99]). But there are two new 
elements as well: i) the specific formulation of the adiabatic condition (which is succinct and 
yet clarifies some subtleties); and ii) the discussion of the regularized Hamiltonian operator 
Hi. 

A. The adiabatic condition 

As explained in Section V B, a choice of a basis of 'generalized positive frequency' solutions 
6^(77) satisfying the normalization condition (5.10) determines a vacuum state, |0), from 
which a Fock space "Hi can be constructed. Since each of these (complex) basis vectors 
efc(f/) satisfies the second order, linear, ordinary differential equation (5.9), any two bases 
efc(f^) and e;i.(?7) are related simply by [93, 94] 

ekiv) = akei,{fj) + Pkelifi) , (6.1) 

where the time-independent Bogoluibov coefficients ak and Pk satisfy the relation — 
= 1. Substituting this equation in the expression (5.8) of the field operator T(a;, r/), 
we find a linear relation between the creation and annihilation operator associated with the 
two families 

A,, = ak A^^ + (31 Al^^ . (6.2) 

This relation shows that, as long as the coefficient are not trivial, i.e. 7^ e*^'= for all 
k, the associated vacua |0) and |0) are distinct. The number of 'under-barred' quanta with 
momentum k that the state |0) contains is given by 

(0|N,HO) := mhfr'A^,^^ = \(3,\' , (6.3) 
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where we have used (6.2) and the commutation relation (5.11) in the last step. The right 
hand side provides the expected number of 'under-barred exitations/particles' with momen- 
tum k in the vacuum |0). Therefore, at the power counting level, it follows that if |/3fcp 
does not fall-off faster than k~^ when k — )■ oo, the total number of 'under-barred quanta' 
in the vacuum |0) would diverge. In this case, |0) would not belong to the under-barred 
Fock space; the two representations of the Heisenberg or Weyl algebra would be unitarily 
inequivalent. 

This inequivalence is related to the large k limit and, as indicated in section VB, can 
arise because so far there is no restriction on the ultraviolet behavior of the basis states 
ek{fj). The physical idea behind the appropriate restriction becomes clearer if one works 
with the variable Xk{v) ■= {(^{v)/'^\^) ek{fj), for which (5.9) becomes 

xl{v)+(k'-~^] XM = 0. (6.4) 



Note that (6.4) reduces to the equation satisfied by the standard basis functions in Minkowski 
space if a" /a = 0. Now, Qa" /d^ is just the scalar curvature of (jab and introduces a physical 
length scale L^fj) into the problem. The form of (6.4) suggests that for modes with large 
momentum, i.e., with k/d ^ curvature would have negligible effect and they would 
evolve almost as if they were in Minkowski space. Therefore, it is natural to impose the 
following regularity condition on the choice of basis functions in the ultraviolet limit: for 
/c/a ^ 1/L, Xkiv) should approach the canonical Minkowski space positive frequency solu- 
tions e*'^*/ \/2k at an appropriate rate. (In the terminology of section VB, we would then be 
restricted to a preferred family of complex structures all of which have the same ultraviolet 
behavior as the canonical complex structure in Minkowski space-time.) This is the crux of 
the idea behind adiabatic condition. 

To make this idea precise, we have to sharpen the required rate of approach. For this, 
one first introduces a set of specific 'generalized WKB' solutions to (6.4) that approach the 
Minkowski space positive frequency modes as /c — )■ oo in a controlled fashion. In a second 
step, one requires that the permissible basis functions Xkiv) — which are exact solutions to 
(6.4) — should approach the WKB solutions to the desired order. 

The generalized WKB solutions Xk^\v) of order are given by [95]: 

\/21^r(^) 

where wl^^\fi) = Wo + Wi + ... + Wn, with 

, 1 d" 2a"a'2 - 2d"^d - 2a5'fi"' + d^d"" 

Wo = k; W2 = - — —] W4 



2k d' " 8W ' 

and Wi = Q if z is odd. (6.6) 



For this prescription to be well-defined, {ff) must be non- negative. For any given smooth 0(7)), this 

can be ensured by going to high enough k. Since it is only the behavior of wl.^\fj) in the /c — > 00 that 

matters for the adiabiticity considerations, if W^!^^ were to become negative in some fc-range, one can 
just suitably modify its form for low k. 
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Note that, because the lower bound on the integral in the phase factor is not fixed, Xk 
well-defined only up to an overall phase (which is time independent but can depend on k). 
Each Xk^\v) is an approximate solution to (6.4) in the sense that, when we operate on it 
by the operator on the left side of (6.4), the result does not vanish but is given by terms of 
the order 0{a/kL^^2)^~^^^^ where the length scale L^v is dictated by d and its //-derivatives 
to order A^. The leading order term of (6.5) corresponds to the positive frequency solution 
in Minkowski space and the rest of the terms are higher order contributions that vanish at 
different rates when {a/kLj^) — >■ 0. Finally, this approximate solution can also be regarded 
as an expansion in the number of derivatives of the scale factor. This is because, since the 
limit refers to (d/kL^), one can either keep L^r fixed and let d/k go to zero as we have done 
so far, or keep d/k fixed and let I/Ltv go to zero, which corresponds to letting the expansion 
rate go to zero, i.e, considering adiabatic expansion. This why, although this method is 
primarily concerned with ultraviolet issues, it is referred to as 'adiabatic regularization'.^^ 

We are now ready to state the 'adiabatic condition' that imposes the desired ultra- 
violet regularity on the basis functions: In the mode expansion (5.8), choose only those 

solutions ek{fj) = 2y^Xfc(^)/^(^) to (5.9) which agree with x\^\v) W terms of order 
{d/kLN)^^^^"^ . More precisely, we require: \xk\ = 1x1^'' I (1 + C){{d/kLN+2))^^'^)-, and 
the same relation should hold if Xk and xi^^ are replaced by their 1st, 2nd, ... (N-l)th 
time derivatives. (Note that the absolute value signs make these conditions well-defined in 
spite of the phase ambiguity in xi^''-) These bases ek{fi) will be referred to as Nth-oidei 
adiabatic solutions, the associated vacuum states will be referred to as the Nth order 
vacua, and states obtained by operating on these vacua by (arbitrarily large but finite) 
sums of products of creation operators as Nth-order adiabatic states. From now on, 
we will restrict ourselves to Nth-order adiabatic states and, for reasons explained in 
section VI B, for most part we will set = 4. There is a large body of literature on the 
notion of adiabatic states and their properties. For further details, see in particular [8, 9, 95]. 

This framework has two important features. 

• The adiabatic condition is only an asymptotic restriction for large kL/d. Therefore, for 
any given iV, there are infinitely many families of solutions ek{fi) which satisfy it. Each 
of these bases defines an adiabatic vacuum and, if two bases are non-trivially related 
(i.e., if the Bogoluibov coefficients are not pure phases for all k), the corresponding 
vacua are distinct. Thus, in striking contrast to the free field theory in Minkowski 
space, there is no preferred vacuum state, nor a canonical notion of 'particles'. Any 
one vacuum appears as an 'excited state with many particles' with respect to another 
vacuum. 

• However, if A^ > 2, all adiabatic vacua — and hence all adiabatic states — lie in the 
same Hilbert space "Hi. This is because if A^ > 2, then ^ < oo, whence any one 
adiabatic vacuum has only a finite number of particles relative to any other. 



Chronologically, the term adiabatic originated from the fact that the first application of this method — 
introduced by Parker in [93] — was to the problem of regularizing the number operator in an expanding 
universe, and the condition that the particle number be an adiabatic invariant was used as a fundamental 
requirement in the construction. 



33 



This completes the specification of the Hilbert "Hi of perturbations we began in section 
VB. While we used the framework originally developed for quantum field theory in classical 
space-times, because we used the dressed effective metric cjab for our background, it follows 
from section IV that "Hi is the Hilbert space of quantum perturbations ip propagating on 
the quantum geometry \I^o- 

Remark: Had we chosen to work with spatial manifolds M with M.^ topology rather than 
T^, Eq (6.3) would be replaced by 

(0|N,HO) := m-%lJm = \M'{27rrS'{0) . (6.7) 

This implies that \/3k\'^ is now the number density of the 'under-barred exitations/particles' 
— with momentum k — in the vacuum |0) per unit volume of M and per unit co-moving 
volume in the momentum space. The 5^(0) (in the momentum space) in (6.7) arises because 
of the infinite spatial volume of M. Thus, in addition to the potential ultraviolet divergence 
that would occur if does not fall-off appropriately for large k, we now also have an 
infrared divergence. Note that this cannot be cured simply by putting an infrared cut-off in 
the k space: the total number of particles created with momenta k within any finite range 
Afc also diverges because of the infinite spatial volume of M. In particular, this divergence 
persists for massive fields as well; it arises because we now have an infinite volume and a 
homogeneous background. This then implies that the Fock representations of the Heisenberg 
or Weyl algebra associated with any two bases efc(f/) and e^(//) are unitarily inequivalent 
unless ak is a pure phase for all k. But physically this infinity is spurious. Therefore, for 
topology, a notion of 'physical equivalence' is more appropriate than that of 'unitary 
equivalence'. To introduce it, recall first that if the topology is one has to introduce a 
fiducial cell C and restrict all integrations to it already at the classical level. (From a physical 
perspective one can choose the cell so that its physical volume is larger than the volume 
of the observed universe.) Two Fock representations of perturbations would be regarded as 
physically equivalent if the vacuum state associated with one contains a finite number of 
'exitations/particles' with respect to the other within any region contained in C. Then, if we 
restrict ourselves to adiabatic vacua of order N > 2, we are ensured of physical equivalence. 
If we require > 4, the expectation values of the regularized stress-energy tensor would be 
well-defined distributions and we can restrict the support of test functions in C. 



B. Regularization of composite operators 

In this subsection we summarize the necessary regularization procedure to obtain physical 
information from the formal expressions of composite operators on the Hilbert space "Hi. 
Adiabatic regularity of the basis modes efc(f/) will provide the necessary control on the 
ultraviolet divergences in the expectation values of composite operators, leading to state 
independent criteria to extract the physical, finite results while respecting the underlying 
covariance of the theory. 

Consider a formal operator 0{x, fj) which is at least quadratic in the field operator and its 
conjugate momentum, and is factor ordered to be self-adjoint. Examples of direct interest 
are the stress-energy tensor and the Hamiltonian. Consider the expectation value of this 
operator with respect to an adiabatic vacuum selected by a basis Xk{v)- Using the mode 
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decomposition (5.8) and the commutation relation (5.11), the vacuum expectation value of 
0{x,fi) can be expressed as a formal sum in the momentum space, of the type 



where there is no x-dependence on the right side because |0) is translationally invariant. 
Generically, the sum would be ultraviolet divergent. In adiabatic regularization the phys- 
ically relevant, finite expression is obtained by subtracting, mode by mode, each term in 
the adiabatic expansion of Ok[xk{v)] that contains at least one ultraviolet divergent piece 
[93, 95] (see also [8, 9]). Thus, if O^^^ is the Nth adiabatic order term in the expansion of the 
summand, if any one part of O'-^^ is divergent, the entire term O^^^ is subtracted (including 
parts that have no ultraviolet divergences). On the other hand, following the criterion of 
minimal subtraction, this procedure is applied only up to that order in the adiabatic expan- 
sion at which the formal expression has divergent pieces. Thus, if the Mth adiabatic order 
term has no divergent part, then nothing is subtracted at order M. 

The most interesting example for us is the stress-energy tensor of gauge invariant tensor 
perturbations since it plays the key role in checking if the truncation scheme is self-consistent, 
i.e., whether or not the back-reaction can be neglected. For classical fields T{x), the expres- 
sion is given by 

Tai, = [VaTVftT- ^^a^r' VeTV.r] . (6.9) 
4/t 2 -' 

At the quantum level the stress-energy tensor is a composite operator of dimension four, 

and in four space-time dimensions we expect ultraviolet divergences up to fourth order in 

the CO- moving momentum k. Let us first consider energy density operator p{x,fj). Given a 

basis ek{fj) = 2v^ Xfc(^)/^(^) of fourth or higher adiabatic order, the formal expression for 

the expectation value of p{x, ff) in the associated adiabatic vacuum is 

(0|p|0)fo™al := -(0|TO|0)formal=^^^X]pfc[Xfc(^)] 

k 

k 

By using the adiabatic expansion of each Xk in the above summand, it is easy to see that 
all the ultraviolet divergences are contained in terms of adiabatic order equal to and smaller 
than four. The zeroth adiabatic order term produces a ^(l//c'^) divergence in co-moving 
momentum k; the second order term a ^(1/fc^) one; and the fourth order term a XKV^) 
one. Therefore, the subtraction terms C^(/c,f/) needed to regularize the energy density are 
obtained from the terms of zeroth, second and fourth adiabatic order in the expansion of 
the summand: 

^, (0) (2) (4) k a'^ 4a'^a" + aa"^ - 2aa' a" 
Thus, the vacuum expectation value of the renormalized energy density is 

(0|P|0)- = ^5-^ E ^P>^\Xkm - C^ik, f,)) . (6.12) 

k 
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Note that the subtraction terms are local in the background geometry and, even more 
importantly, are state independent. The expectation value in any fourth order state in "Hi 
is computed by the same procedure, using the same subtraction terms. 

The only other independent degree of freedom in the stress-energy tensor in a homogenous 
and isotropic background is the trace T (or the pressure p = (p + T)/3). The corresponding 
vacuum expectation value can be calculated using results given in [99]: 

k 

where the 4th adiabatic order subtraction terms are 

^rj. —ol'^ + aa" Qa''^a'' — ?>aa'''^ — Aaa'a' + a?a''' , ^. 

^ ^) = 8W • ^^-^'^ 

Note that both vacuum expectation values, (6.12) and (6.13), are constant functions on M 
as they must be given that the vacuum is translationally invariant. jf the vacuum were 
replaced by a generic 4th order adiabatic state, this constancy will not hold. In that case 
the expectation value would be distributions in x. Finally, by themselves these expectation 
values only provide a quadratic form on the Hilbert space "Hi. However, recent results show 
that they are the expectation values of an operator valued distribution Tab on "Hi [96]. 

A natural question now arises: Does the vacuum expectation value of the renormalized 
stress-energy tensor satisfy covariant conservation Va(T"'')rcn = with respect to the dressed 
effective background cjab'^ In any homogenous and isotropic background, the conservation 
equation reduces to + |^(4(/))rcn+ (T)rcn) = 0. The formal, unrenormalized expressions 
for (p) and (T) satisfy this relation mode by mode as a consequence of the wave equation 
satisfied by the Xk- The adiabatic subtraction terms also satisfy the conservation equation 
mode by mode. Thus, one can directly verify that the expectation value of the renormalized 
stress-energy is indeed conserved (see, for instance, [97]). 

This adiabatic regularization of Tab also has the desired properties enunciated in Wald's 
axioms [92]: it reduces to the standard normal ordering in the fiat space-time limit; the 
subtractions terms are constructed from local information in the background geometry, and, 
as already noted, the renormalized stress energy is conserved; Va(T'^^)ren = 0. 

Remark: As noted in section III, physically it is appropriate to introduce an infrared 
cutoff by absorbing into the background the modes whose physical wave length is larger than 
the physical radius of the observable universe. Since physical wave lengths scale linearly 
with the scale factor a, this can be achieved in a time independent fashion by imposing 
a cut-off, kin, in the co-moving wave number k. As is clear from the above discussion 
of conservation of stress-energy, the vacuum expectation value of the new renormalized 
stress-energy tensor will also be covariantly conserved. Furthermore, by construction, the 
renormalized energy density and pressure will again be constant on M. 

We will conclude this section with renormalization of the Hamiltonian operator Hi that 
generates the dynamics of perturbations in conformal time fj. The form (3.31) of the classical 



As a consequence, while the summands on the right side of these expressions can be interpreted as the 
contribution of the mode k to (p(a;))ren and (r(a;))ron, they are not the Fourier transforms (/5fc)ion and 
{Tk)ren of thcse quantities in the x-space. 
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Hamiltonian and the fact that the lapse corresponding to the conformal time is 
that the Hamihonian operator has the following formal expression: 

TT X ^ 4 K . ^ 12 a J 2 \ tn" 1 2 ~ 4 f j3 ^ 

-ni formal = —o; / , "ry iPul + k Tr = a d X Pformal • 

k 

Therefore, the renormalized Hamiltonian is given by 

i^l,ren = i | + ^ fc^ 1 - /iV^ ^''(fc, r)) . (6.16) 

k 

Our discussion of energy density immediately implies that the expectation values of Hi^rcn 
are well defined on any 4th order adiabatic state. Furthermore, by taking its commutators 
with T^ipj: we recover the Heisenberg equations of motion (4.18) satisfied by 7^. However, 
we do not have a proof that Hi^^cn is a self-adjoint operator (in the precise sense that it 
is densely defined and its domain equals that of its adjoint. In particular, the strategy of 
using the Friedrich extension does not work because the operator is not positive definite.) 
Although we expect this to be the case, we also note that, quite generally in cosmology, 
proving self-adjointness involves non-trivial subtleties and technicalities especially because 
the Hamiltonians have a non-trivial time dependence. 

Remark: In the case when M is topologically M^, one continues to use the mode by 
mode subtraction strategy but one has to replace (1/-^^)^^ by l/(27r)^ J d?k. Again, the 
renormalized stress-energy is conserved, the subtraction terms are local and the prescription 
reduces to the standard normal ordering in Minkowski space-time. Again it is physically 
appropriate to introduce an infrared cut off in the co-moving /c-space and conservation 
of stress-energy persists after imposing this cut-off. In addition, as shown in Ref. [98, 
99], in FRLW space-times this adiabatic regularization is equivalent to the point-splitting 
Hadamard renormalization. Therefore although the procedure is not manifestly covariant 
because of the mode by mode subtraction, the result is fully covariant w.r.t. gab- From 
a fundamental quantum consideration, cjab is only a convenient mathematical construction 
and the true quantum geometry is encoded in ^ o- Still, it is desirable that the effective 
description have space-time covariance. 



C. Narrowing down initial conditions 

As discussed in section VIA, our quantum Hilbert space T-ii of perturbations does not 
admit a preferred vacuum state. In quantum field theories on Minkowski or de Sitter space- 
time, underlying isometrics serve as powerful tools to single out preferred states. Our quan- 
tum FLRW background o is invariant under the 3-dimensional translation group acting on 
M = T^. It is then natural to seek preferred states in T-ii by demanding that they also be 



Counter terms now appear under an integral sign. Results of Ref. [98, 99] imply that, one can use 
this integral form of counter terms also in the spatially compact cases. This procedure is justified on 
the grounds that the ultraviolet regularization should not be sensitive to the global topology and has, 
furthermore, some advantages. 
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translationally invariant. As discussed in section V, this condition leads us to the infinite 
dimensional space of vacua which, in view of our discussion in VI B, we now require to be of 
4th adiabatic order. If one requires the state to satisfy these properties initially, i.e. at the 
bounce, then they are satisfied for all times. This is the family of 'preferred' states selected 
by the symmetry and regularity requirements. As we will discuss in detail elsewhere, this 
choice is well-suited to formulate a quantum version of Penrose's Weyl curvature hypothesis 
[100]. Furthermore, in the infiationary context, one can motivate this choice using phys- 
ical considerations based on properties of the new repulsive force with origin in quantum 
geometry that dominates the dynamics at and near the bounce [1, 2]. 

At a fundamental level, we allow any of these 'vacua' as initial conditions at the bounce 
provided, of course, the energy density (p) in the perturbations at the bounce is negligible 
compared to the energy density (po) ~ 0.41ppi in the background quantum geometry "^g. 
However, for specific calculations and especially detailed numerical simulations, one has to 
work with specific states. Are there then especially convenient vacua to work with? The 
adiabatic procedure summarized in section VIA provides a strategy to select an 'obvious' 
candidate, provided one fixes an instant of time fjo. Recall that to ensure 4th order adia- 
baticity we required that, for kL^/a S> 1, the basis functions Xk must agree with a specific 
approximate solution xC' defined in (6.5) at least up to terms of adiabatic order four. There- 
fore, given an instant fjQ of time, we can construct a 'natural' basis xt'^ by asking that it 

has the same initial data at that time as x^k' ■ Thus the idea is to ask for solutions X^k'^iji) 
to the exact evolution equation (6.4) which satisfy 

(^o) = xt\m) ; and d^xt\fl,) = dr-,xf{m) ■ (6-17) 

Since x^k^ ^'^^ only approximate solutions to (6.4), they will not agree with xt^^ at any other 
time fj ^ fjQ. Nonetheless, X°k''{v) are automatically of 4th adiabatic order for all times 
fj. We will call the associated vacuum state the obvious 4th-order adiabatic vacuum at the 
time fjQ. In LQC, the preferred instant of time required in this strategy is provided by the 
bounce. 

Note, however, that for this strategy to work, the quantity W^ifjo) in (6.5) must be 
non- negative since it appears under a square-root in the expression of X^k\^o)- Now, as we 
pointed out in III, there is a natural physical infrared cut-off kj^ provided by the radius of 
the observable universe. Typically, for k > /cir,, i.e, for the relevant k, wjf^ is positive. But 
if km is too low for this to happen, as explained in section VI A, one has to suitably modify 
the form of W^lfjo) for low k. In this case, even for a fixed value r/o of time, there would 
be an ambiguity in the choice of the basis Xk for low k, and hence also in the resulting 4th 
adiabatic order vacuum jO)^!;"^. 

This vacuum state |0)^|^^ is especially convenient to work with in numerical simulations. 
Therefore it serves as a technically powerful tool to establish the viability of conjectures — 
e.g., that, in the infiationary context, there exist quantum states for which the back-reaction 
can be ignored — and to probe qualitative features of quantum dynamics. However, even 
when there is a preferred instant of time, such as the bounce time in LQC, and the infrared 
cut-off is sufficiently large for the state to be unique, there is no physical principle that 
singles out this vacuum over other 4th adiabatic order vacua. A better strategy would be 
to narrow down the choice of vacua by adding external inputs suggested by the physics of 
the problem. In the present analysis there is a natural avenue along these lines. Since the 
test field approximation plays a central role, one could ask: Are there 4th adiabatic order 
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vacua for which the expectation value of the renormahzed energy density is exactly zero at 
a given time? The answer to this question is in the affirmative and, although this condition 
does not single out a unique state, it reduces the choices considerably. This issue will be 
discussed in detail in a separate publication [101]. 

D. A Criterion for self-consistency of the truncation scheme 

Our results in sections V - VI C provide us with a specific quantum theory, corresponding 
to the truncated phase space FTrun. This self contained mathematical framework enables one 
to describe quantum dynamics of the truncated system, once the initial state is specified, 
say at the bounce. The key question now is whether the truncated theory has an interesting 
domain of validity. The basic assumption behind truncation is that the back-reaction due 
to the stress-energy tensor Tab of the first order fields can be neglected (compared to the 
background stress-energy tensor T°^) during the dynamical phase of interest. Already in 
the classical theory, it is clear that generic first order perturbations violate this condition. 
The interesting question there is rather the following: Is there a sufficiently large subspace 
in FTrun on which this condition is satisfied? More precisely, if we restrict ourselves to an 
instant of time and choose perturbations which do meet this condition, does the condition 
continue to be satisfied over the period of evolution of interest? We will now formulate the 
analogs of these questions in the quantum theory. 

Somewhat surprisingly, in the quantum theory one encounters some novel features. First, 
in the classical theory the only 'preferred' time to specify the initial conditions would be 
the big bang which is singular. In the quantum theory, by contrast, the bounce provides 
a natural time r^^ for this task. The second feature is equally significant but also more 
subtle. Because of the underlying symmetries of the background, it is natural to require 
the initial state ip at time r^^ to be translationally invariant. An immediate consequence is 
that the expectation value (Tab) of the renormahzed stress-energy tensor operator is then 
homogeneous for all times}^ This property has the important consequence that the second 
order perturbations are again homogeneous and isotropic since their sole source is {Tab)- 
Thus, thanks to the symmetry of our initial state if) at the bounce, the back reaction can 
change only the zeroth order, homogeneous fields, changing the total state \E'o®V' to a nearby 
state of the type {^o + ^^o) ® V"- Our initial conditions on if) at fj = fj-^ guarantee that the 
shift 5'^ o at this initial time is negligible. The key question then is whether it continues to 
remain negligible under time evolution. 

On general grounds one would say that the answer is dictated by the time dependence 
of (Tafe). But in the detailed mathematical framework, it is the second order truncation §2 
of the Hamiltonian constraint that determines the change 6^ o in the background quantum 
geometry and in this equation it is only (p) — rather than the full {Tab) — that enters as the 
source. Let us suppose that the energy density (p) is negligible compared to the background 
energy density (po) from r^^, until some time fjQ of interest. Then it follows that 5^o would 
continue to be negligible from fj-Q to fjQ. It may seem somewhat surprising at first that one 
does not have to require explicitly that the other independent component T of Tab should 



Note that this property holds even though the perturbations operators 77) themselves are purely 
inhomogeneous\ In the classical theory, there are no perturbations T{x, 77) which are purely inhomogeneous 
but their energy density is purely homogeneous. 
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be small. This is because we are requiring that (p) be negligible compared to (po) not just 
initially but for all times between fj^ and fjo . The difficult part of the calculation is to check 
that the evolution of ip is such that this condition does holds. If it does, then ® would 
provide a solution in which the back-reaction is negligible even in the Planck era. This 
would then be a self-consistent solution to the quantum truncated theory. Thus, a sufficient 
condition for self- consistency of the truncation approximation is that the energy tensity in 
the quantum perturbations should remain small compared to that in the background from fj-^ 
to fjo. 

In [2] we will use detailed numerical simulations to show that the inflationary scenario 
does admit states \l/o ® in which (p) is negligible compared to the background (p^) for 
all times between the bounce and the onset of the slow roll inflation. Furthermore, given a 
state satisfying this condition, we will show that there is an open neighborhood of tp such 
that the condition continues to be satisfied by \l/o ® ^Z' for all ijj in this neighborhood. Thus, 
there is a rich class of states that provide self-consistent solutions to the truncated quantum 
theory, demonstrating that the standard inflationary scenario admits a consistent extension 
all the way back to the big bounce. 

But the general framework constructed in this section is not tied to inflation. It provides 
the technical machinery that is needed to check if any given paradigm, based on general 
relativity and ffist order cosmological perturbations, admits a self- consistent extension to 
the Planck regime. More precisely, it would enable one to address the following question: 
In this paradigm, does the quantum theory admit solutions in which the back reaction can 
be neglected throughout the period of interest, including the Planck era? 

VII. SUMMARY AND DISCUSSION 

In the last four sections, we developed an extension of the standard cosmological per- 
turbation theory to include the Planck regime of LQG. The strategy was to ffist truncate 
classical general relativity coupled to a scalar field to the sector commonly used in the cos- 
mology of the early universe — FLRW space-times and linear inhomogeneous perturbations 
thereon — and then construct the quantum theory of just this sector using LQG techniques. 

Already in the truncation of the classical Hamiltonian theory there is a subtlety that 
has often been overlooked in the LQC literature: while the dynamics on the homogeneous 
sector is generated by a (Hamiltonian) constraint, that on the full truncated phase space 
TTrun is uot. Indeed, because the dynamical vector field X"' on FTmn fails to Lie drag the 
full symplectic structure ^Trun = fio + ^i on FTrun, it is not generated by any Hamiltonian. 
Rather, Fttuii is the normal bundle over Fq — where the base space F,, is regarded as the 
homogeneous, isotropic subspace of the full phase space F of general relativity — and X" is 
the lift to Fxmn of the Hamiltonian vector field on F^, induced by the full Hamiltonian vector 
field on F. In the classical theory, this subtlety can be ignored if one works with the space of 
solutions (rather than the phase space), as is common in the standard cosmology literature. 
But it becomes important for passage to quantum theory if one wishes to treat both the 
perturbations and the background quantum mechanically. Then there is no conceptual 
justification for trying to construct dynamics for the full truncated system by imposing a 
quantum constraint. One can do this only on the homogeneous sector, and one then has to 
'lift' this quantum dynamics to the full Hilbert space just as in the classical theory. 

Having constructed the dynamics of gauge invariant variables on the truncated phase 
space, we then used LQG techniques to construct quantum kinematics: the Hilbert space 
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"Ho of states of background quantum geometry, the Hilbert space "Hi of gauge invariant 
quantum fields Q, T representing perturbations and physically interesting operators on both 
these Hilbert spaces. The imposition of the quantum constraint on the homogeneous sector 
leads one to interpret the background scalar field as a relational or emergent time variable 
with respect to which physical degrees of freedom evolve. Furthermore, the background ge- 
ometry is now represented by a wave function \l/o which encodes the probability amplitude 
for various FLRW geometries to occur. The physically interesting wave functions \I^o are 
sharply peaked, but the peak follows a bouncing trajectory, not a classical FLRW solution 
that originates at the big bang. In addition, has fluctuations about this bouncing tra- 
jectory. Quantum fields Q, T, representing inhomogeneous scalar and tensor perturbations, 
propagate on this quantum geometry and are therefore sensitive not only to the major de- 
parture from the classical FLRW solutions in the Planck regime, but also to the quantum 
fiuctuations around the bouncing trajectory, encoded in \l/o. Therefore at first the problem 
appears to be very complicated. However, a key simplification made it tractable: Within 
the test field approximation inherent to the truncation strategy, the propagation of Q, 1~ 
on the quantum geometry \l/o is completely equivalent to that of their propagation on a 
specific, quantum corrected FLRW metric gab- Although h does appear in its coefficients, 
this 'dressed, effective metric' gab is smooth and allows us to translate the evolution of Q, T 
with respect to the relational time to that in terms of the conformal (or proper) time of gab- 
Furthermore, away from the Planck regime, gab satisfies Einstein's equations to an excellent 
approximation. In this sense, the standard quantum field theory of Q, T emerges from the 
more fundamental description of these fields evolving on the quantum geometry ^ „ with 
respect to the relational time 0. This exact relation between quantum fields Q, T on the 
quantum geometry and those on the dressed, effective geometry of gab enabled us to carry 
over adiabatic regularization techniques from quantum field theory in curved space-times to 
those on quantum geometries ^o- Together, all this structure provides us with a well-defined 
quantum theory of the truncated phase space we began with. 

This framework has a broad range of applicability because scenarios of the early universe 
are often based on linear perturbations on FLRW backgrounds. Our construction provides 
an avenue to extend them all the way to the quantum gravity era because quantum pertur- 
bations now propagate on a quantum geometry which is completely regular, with specific 
upper bounds for curvature and density in the background. We can therefore use the new 
framework to re-examine the 'trans-Planckian issues' encountered in these scenarios. Note 
first that, the truncated theory under consideration here allows modes with trans-Planckian 
frequencies. There is no obstruction because the quantum geometry underlying LQG is sub- 
tle: In particular, while there is a minimum non-zero eigenvalue of the area operator, there is 
no such minimum for the volume or length operators even though their eigenvalues are also 
discrete. The real danger is not the existence of such modes but rather that the energy density 
in these modes may not be negligible compared to that in the quantum background geometry. 
If this occurs, our quantum theory of the truncated sector would not be viable. Whether this 
can happen is a very non-trivial issue especially in the Planck regime immediately following 
the bounce. Heuristically, if the state has just a few excitations each carrying say, 10^ times 
the Planck energy in a cm^ volume, there would be no difficulty (since the energy density 
would be negligible). If on the other hand there is one such excitation per Planck volume, 
our truncation approximation will fail. Then we cannot neglect back-reaction. This is not 
an impasse to quantum theory as such, but the proper treatment of such states will have to 
await full LQG. 
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The key question then is whether the test field approximation underlying this truncation 
scheme is satisfied. A priori this is a difficult issue and, to our knowledge, had not been 
considered in the literature because even to formulate this question precisely one needs the 
notion of the renormalized stress-energy tensor on the Hilbert space "Hi of quantum pertur- 
bations propagating on the quantum geometry of \E'o- In our framework this was provided 
by 'lifting' the adiabatic techniques of [95, 99] to quantum fields on quantum geometries "^o- 
Specifically, by appealing to symmetry principles and regularity requirements, we argued 
that it was appropriate to focus attention on those states ^'o'^V^ G 'Ho<S)'Hi of the combined 
system for which: 

1) is invariant under the translational symmetry of "^o] 

2) ip is a 4th order adiabatic state w.r.t. cjab] and, 

3) at the bounce, the energy density (p) in the state ip is negligible 
compared to the energy density (po) in the background. 

We then showed that the truncation approximation is self-consistent if (p) continues to 
remain negligible compared to (po) from the bounce time f/^ to a late time fjo of physical in- 
terest (e.g., when radiation decoupled from matter). In particular, our argument shows that 
the full stress-energy tensor is not needed; this significantly simplifies the task of performing 
numerical simulations that are needed to check self-consistency. 

As noted in section I, this criterion does not imply that truncated solutions are necessarily 
close to exact solutions because the sum of all higher order effects need not be negligible. 
However, in practice such criteria are generally regarded as sufficient for truncations to be 
trustworthy. Indeed, this philosophy governs the entire theory of cosmological, stellar and 
black hole perturbations in general relativity as well as perturbative calculations in quantum 
ffeld theory. In the same spirit, our self consistency criterion can be used to test viability of 
the ffrst order truncation in the studies of the very early universe that include the Planck 
regime. 

In the next paper [2] we will use this criterion in the context of inffation. We ffrst extend 
the general framework of this paper slightly to incorporate the (l/2)m^0^ potential and 
spatial topology. (The value of m is fixed by using the 7 year WMAP data [102, 103].) We 
motivate the initial conditions — called 'quantum homogeneity' at the bounce — and carry 
out detailed numerical simulations using the 'obvious' 4th order adiabatic vacuum for the 
initial quantum state ip of perturbations. They show that the back-reaction of perturbations 
remains negligible over the 11 orders of magnitude in matter density and curvature, from 
the big bounce until the onset of slow roll. Furthermore, the power spectrum at the end of 
inffation turns out to be very close to that obtained in the standard inffationary scenario 
and is thus compatible with the WMAP observations. By varying initial conditions for 
the background (within computational feasibility) we show that these results are robust. 
Furthermore, we show that self-consistency is preserved if the initial state is chosen to 
be in a neighborhood of the 'obvious' 4th adiabatic order vacuum. Taken together these 
results establish existence of self-consistent extensions of the inffationary scenario to 
the Planck regime. Finally, there is a small range for the value 0b of the background 
inffaton ffeld at the bounce for which the quantum state at the onset of inffation differs 
sufficiently from the Bunch Davies vacuum assumed in standard infiation to give rise to 
non-Gaussianities that could be measured in future observations along the lines of [104-107]. 

We will conclude by pointing out a direction for significant improvements and extensions 
of this framework. We began with a truncation of general relativity coupled with matter. 



42 



that is well-suited for cosmology of the early universe. In the passage to the quantum the- 
ory, for the homogeneous sector we used LQC framework, rooted in the quantum geometry 
underlying LQG. This was crucial for the resolution of the big bang singularity and the 
subsequent quantum dynamics in the Planck era. On the other hand, to facilitate com- 
parisons with the standard cosmological literature, we used a Fock-type representation for 
perturbations Q, T. As in the Gowdy models [25-29], this is a well-defined and internally 
consistent quantization. But it would be more satisfactory to use a 'polymer-type' represen- 
tation rooted in LQG also for perturbations not only for aesthetic reasons but also because 
it would provide sharper guidelines to relate the truncated quantum theory to full LQG. 
Therefore, let us first ask: Would this change of representation make a qualitative difference 
in the results? The following heuristics lead us to believe that the answer is in the nega- 
tive. Note that Q, T represent perturbations and, in any self-consistent solution ® tp, the 
energy density in the perturbations is negligible compared to that in the background}^ On 
general grounds one would expect that, in a viable LQG representation of states capturing 
this physics, their dynamics and properties would be well approximated by our G "Hi. 
As a concrete illustration, we can use the following simplistic strategy to aid intuition. Use 
the background structure available in the truncated sector to decompose the perturbations 
Q{x), T{x) into Fourier modes thereby representing these fields as an assembly of harmonic 
oscillators, and then imagine using the polymer representation of these oscillators to con- 
struct the LQG Hilbert space for these fields. Then one knows [108, 109] that for low 
energy states the results of the polymer and the standard (Fock-type) quantization are in 
excellent agreement. This suggests the 'hybrid' approach is viable from practical or phe- 
nomenological perspective. Moreover, it is well suited to bridge quantum field theory on 
quantum geometries with the well-established quantum field theory in curved space-times. 

However, from a fundamental perspective it is highly desirable to systematically extend 
this framework by replacing T-ii with an appropriate Hilbert space that descends from LQG. 
In particular, such an extension will enable one to arrive at the regularization and renor- 
malization procedure 'starting from above' i.e., from full LQG considerations. By contrast, 
in this paper, we have introduced this procedure 'starting from below', i.e., from quantum 
field theory in curved space-times. Put differently, our primary goal of this paper is to 
carve out a path to extend the cosmological perturbation theory to the Planck era. The 
emphasis has been on showing that there does exist such a framework with a number of 
desirable mathematical properties and physical features which, moreover, is well suited for 
phenomenological applications. But from a fundamental LQG perspective, it can and should 
be related to LQG even more closely. 
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Appendix A: Truncated dynamics: A simple example 

In this appendix we will illustrate the truncation procedure of section III using A$^- 
theory. This example is simple enough to perform explicit calculations that bring out the 
main conceptual subtleties — in particular the differences between dynamics of the exact 
and truncated theories — which are sometimes overlooked in the cosmology literature. 

1. Space-time framework 

As in the main text we assume that the space-time M is topologically M x M, where 
the Cauchy surfaces M are topologically 3-toruses T^. But the space-time metric rjab is now 
assumed to be flat (with signature -,+,+,+). Denote by S the space of suitably regular 
solutions to 

□$_^2^_A$3 = o (Al) 

and by So its subspace consisting of spatially homogeneous solutions. We are interested in 
a small neighborhood of So in S. In this neighborhood, it is convenient to consider curves 

2 n 

$ [e] (f , t) = <P{t) + e^^'^ + ^^'^ {x,t) + ... + -^^-\x,t) + ... (A2) 

parameterized by e g] — 1, 1[, which pass through So at e = 0. The (x, t) are to be 
thought of as the nth order, inhomogeneous perturbations on the homogeneous solution 
4>{t). Since we are interested in curves that move away from So, to avoid redundancy, 
without any loss of generality we will assume that the first order perturbation ip^^\x,t) are 
purely inhomogeneous, i.e. that 

Jdvip^^\x,t) = \/t. (A3) 

Here in what follows all the integrals are over M and dv is the natural volume element 
thereon. Substituting the expansion (A2) in (Al) and matching coefficients of e" for each n 
we obtain a hierarchy of equations: 

+ + = 0, (□ - /i^ - 3A02)v?(^) = 0, 

(□ - /i^ - 3A02)y,(2) = 6A0(y.(l))^ (□ - /i' - 3A02) ^(3) = 6A ((^(i))^ + 3^^^'^^^^^), 

(A4) 

Note that the first equation on is non-linear but an ordinary differential equation (ODE) 
and each of the subsequent equation on (yj^"^ is a linear partial differential equation (PDE) in 
with sources containing lower order fields, already determined by solving the previous 
equations in the hierarchy. Thus, the task of solving the non-linear PDE (Al) is reduced 
to solving one non-linear ODE and a succession of linear PDEs.^^ The idea is that an 
approximate solution to the full problem can be obtained by truncating the series to the 



Note that the equation on cp^^-* has no source terms. Therefore, if the initial data is purely inhomogeneous, 
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appropriate order, i.e., by ignoring terms of 0{e^~^^), say. As usual, while the meaning of the 
approximation in terms of the smallness parameter e (which could be tied to the coupling 
constant A for physical reasons) is clear, the truncated series can be a good approximation 
to the full solution $ only if </?*^"'^ remain small compared to 0. It is important to note that 
this scheme of obtaining approximate solutions is distinct from an alternative procedure that 
appears to be used often in the cosmology literature (although sometimes only implicitly). 
That strategy corresponds to defining a field 5$ via $(x, t) = (f){t) + 5$ and solving the 
equation 

□5$ - fi'^6<^ - (3A0^)(5<I> - (3A0) (5$^ - X6<l>^ = (A5) 

To the linear order this strategy agrees with (A4), but not to higher. Indeed, already at 
the second order, one now has to solve a non-linear PDE, with a function as a coefficient, 
which is in some ways more complicated than solving the original (Al). 

2. The Hamiltonian framework 

For the A<l>^ system, one could pass to the quantum theory directly from the classical 
space-time formulation sketched above. However, general relativity is a background inde- 
pendent theory and the generalized Dirac quantization strategy followed in LQC requires 
us to pass through the Hamiltonian framework. Therefore we will now illustrate how the 
truncation procedure of section A 1 works in the phase space language. 

The full phase space F is spanned by pairs ($(x),H(a;)) on the 3-manifold M which is 
topologically T'^ and its homogeneous subspace Fq is spanned by real numbers (0,p(</,)). The 
symplectic structure is given by: 

n{Si, 62) = Jdv 62II - 62^ 5iH] , (A6) 

where 6 = (5$, 611) denotes tangent vectors to F. The corresponding Poisson brackets 
are also the familiar ones: {$(xi), H(a;2)} = S{xi,X2)- Dynamics is generated by the 
Hamiltonian H: 

iJ($, U) = ^Jdv [H2 + + /i2$2 + A $4] _ (^^7) 

Note that although F is infinite dimensional just as the solution space S, and Fq is 2- 
dimensional just as So, there is a key difference: F and Fq are vector spaces unlike S and 
So- As we will see the ability to define constant vector fields on F plays an important role in 
defining truncated dynamics. (This structure is also available on the phase space of general 
relativity and used in the main text to define truncated dynamics.) 

Again, we are interested only in a small neighborhood of the homogeneous subspace 
Fo of F. Therefore we are led to consider 1-parameter family of curves, parameterized by 



so is the solution. Because equations of (/^'"^ for n > 1 have source terms which are non-hnear in the 
lower order fields, we cannot demand that they be purely inhomogeneous. Note also that, for n > 1 there 
is freedom in adding a solution to the homogeneous equation. This can be removed, e.g., by choosing 
retarded solutions to the inhomogeneous equations. 
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$[e](f) = <i) + eLp^'^\x) + . . . + — v?^"^^) + • • • and 

n[e](f) = ^ + eTi^^\x) + ...+ ^7r(")(f) + . . . (A8) 

where Lp^^^x)^^^^^^) are purely inhomogeneous, and Vo is the volume of the 3-torus M. 
By truncating the series at nth order we obtain the phase space that is appropriate for 
describing the homogeneous solutions together with first, second, ... nth order 

perturbations propagating thereon. Since in the main text we focus on just the first order 
perturbations, let us do the same here. 

Thus, let the truncated phase-space Fxrun consist of a doublet of pairs of canonically 
conjugate fields lp^^\x),'k'^^\x)) where 4>,P((t,) are real numbers (representing ho- 

mogeneous fields on M) and ip^^\x),TT^^\x) are purely inhomogeneous fields on M. Thus, 
TTrun = To X Fi. The symplectic structure on FTrun is given by fixrun = fio + ^i, where 

^o{h.52) = 5i052Pw - 5205iPw ^i{5u52) = Jdv [5i^(i)527r(i) - 62ip'^'^6i7rW] (A9) 

which yields the Poisson brackets 

{</., = 1, {^^'\xi), 7r^'\x2)} = [S{xi,X2) -^] = S{Xi,X2) (AlO) 

where S{xi,X2) is the Dirac delta-distribution restricted to the purely inhomogeneous fields 
on M.^^ By decomposing fields $(x),n(x) in F into purely homogeneous and purely inho- 
mogeneous parts, one can readily see that the full pase space (F, Q) is naturally isomorphic 
to the truncated phase space (FTrunj ^Tmn)- However, the physical meaning of the inho- 
mogeneous fields is different in the two cases and, more importantly, the dynamics is very 
different. 

Geometrically, FTrun is the normal bundle over Fq (since purely inhomogeneous fields are 
orthogonal to the purely homogeneous ones in the norm of the space of functions on M.) 
To obtain the dynamical flow on Fxrun, let us begin with the full, non- linear Hamiltonian vec- 
tor field Xh on a small neighborhood of Fo in F. Since the exact Hamiltonian flow generated 
by Xh is tangential to Fq, the equations of motion on Fq are just the restrictions of those 
on full F to the homogeneous sector: = P(^)/Vo, P[ci,) = —{fi'^(j) + \(l)^)Vo. To obtain the 
equations of motion on the (inhomogeneous) tangent vectors {lp^^\ ii^^^) , we first note that 
the Hamiltonian flow Xh on F naturally drags these tangent vectors along dynamical tra- 
jectories on Fq. To obtain the 'dot', however, we need to compare the image {(p^^\Tr^^^)\t+st 
at {(l),p(^^))\t+st with the original tangent vector {lp^^\ n^^'')\t at the point {(j),p[^))\t of Fo. 
This can be trivially accomplished because F has a vector space structure. The resulting 
equations of motion are: ^^^\x,t) = 7i^^\x,t), 7r^^\x,t) = [D^ - fi^ -3X^'^{t)]<p^^^ {x,t). The 



The term 1/Vo is necessary simply because the fields ip^^'){x) and t:^^'>{x) are purely inhomogeneous. For 
example, if we integrate the left side of the Poisson bracket between ip^^^ and tt^^^ over xi (or X2), we get 
zero and 1/Vo terms assures that the right side also vanishes. 



46 



resulting dynamical vector field Xoyn on Fxrun is given by: 

^Dyn := (0,P(</,); 0^^\7r(^) 



(All) 



Thus, the full dynamics on F induces a well-defined flow X^yn on FTrun- Furthermore, 
this dynamical vector fleld can be expressed in a form adapted to symplectic geometry. 
Using Greek letters to denote the abstract indices labeling tangent vectors to F, X^^^ can 
be expressed as 

^Syn = K^dpH, + 9.fdpH^ (A12) 

where 



i7i(y,W,vr«) = ^/rft;[(7r«)2+ +/x2(y,W)2 + 3A(02)(y,(i))2]. (A13) 

Note however that this is not a Hamiltonian flow on Fxrun because it is not of the form 
[nf + nf)dpH^,^^ for any function -ffxrun on Fxrun- The obvious candidate, i^Trun := 
Ho + Hi does not work because Qf^^d^Hi ^ (since Hi depends not only on (v?^^-*, vr*^^)) but 
also on (f))}^ More generally, one can verify that the fact that Hi depends on <^ implies that 
the Lie derivative of l^Trun := + f2i by the dynamical vector fleld Xoyn on Fxrun does not 
vanish. 

In the classical theory, the fact that we have a well-deflned dynamical flow on Fxrun 
suffices. However, the fact that the flow is not Hamiltonian introduces new features in 
the transition to quantum theory. Keeping the quantum perspective in mind, one can 
rephrase the classical djTiamics as follows. Since the homogeneous solution is to be regarded 
as the background and inhomogeneities as perturbations, we can first restrict ourselves 
to the homogeneous part of the phase space (Fq, VLo) and note that, on it, the dynamics is 
indeed governed by a true Hamiltonian fiow, generated by Ho- Fix any dynamical trajectory 
0(t),P((^)(t) in To- To specify how perturbations propagate on this background solution, 
we need to lift this trajectory to the normal bundle, Fxrun- This is precisely what the 
remaining part, Vt^^dp Hi, of the dynamical vector field X^^^ does. Given any tangent vector 

{ip^\ TTo"^'') at a point, say {(f){to),P((i)){to)), along the given dynamical trajectory, the orbit 
of Xgyjj through the point {(f){to),P((f,){to)', fc^\to),'n'c^\to)) of FTrun specifies the dynamics 
of the perturbation tt^^^) on the background trajectory (0(t),p(0)(t)). The fact that 

we have lifts to Fxrun of orbits in Fq trivially implies that, even though different choices of 
the initial {lp^^\7i^^^) define distinct orbits (describing dynamics of various perturbations), 
all these orbits in Fxrun project down to the same orbit on F^. This is just a refiection of 



Had we simply defined dynamics using i/xrun ■— He 
of the baek reaction of the first order perturbation 
because this dynamics ignores the inhomogeneous 
which is of the same order in the e expansion. 



+ Hi , it would have included the homogeneous part 
(i^(^',7r'^)). This would not be physically consistent 
part of the back-reaction of the same perturbation 
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the fact that the dynamics defined by X^^^ neglects the back-reaction of the perturbation 
on the homogeneous background. 

Remark: To obtain the truncated dynamics ^Trun we used the vector space structure of 
r to transport the vector (v^^^-*, 7r(-^))|t+5t at (0,P(</,))|t+5t to the point {(j),p(^^))\t. However, 
one can carry out this comparison more generally, e.g., if there is a natural fiat connection 
to transport vectors from one point of the phase space to another. This is the case when 
r is a cotangent bundle over an affine space (as in LQG) or over a convex subset of a 
vector space (as in the ADM framework of general relativity). Therefore the procedure of 
inducing dynamics on the truncated phase space from the exact dynamics in a neighborhood 
of the homogeneous sector of the phase space goes through also in the cosmological context 
analyzed in the main text. 
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